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INTRODUCTION 

This  report  documents  changes  and  additions  to  a  previous  microwave  calibra¬ 
tion  system  report  (Ref.  1).  In  brief,  the  changes  made  to  the  old  system 
include  the  system  controller  and  software  drivers  for  the  devices  on  the 
GPIB  (IEEE -488)  bus,  and  the  additions  are  the  data  reduction  software  used 
to  generate  tables  and  graphs  from  the  calibration  data. 

This  report  first  discusses  the  reason  for  computerizing  the  calibration  pro¬ 
cess.  Second,  the  hardware  components  necessary  to  complete  the  calibrations 
are  listed.  Finally  the  programs  used  to  obtain  and  reduce  the  raw  data  are 
discussed. 


PURPOSE 

The  primary  reason  for  using  a  computerized  calibration  system  was  to  reduce 
the  amount  of  time  required  to  obtain  the  calibration  and  get  useful  infor¬ 
mation  from  it.  Initially,  the  frequency  response  was  obtained  from  an 
oscilloscope  and  the  data  read  from  a  Polaroid  photo.  The  process  of 
obtaining  the  photo  was  quick,  but  reading  the  data  was  time  consuming  and 
many  errors  could  be  introduced.  When  the  time  for  each  calibration  was  then 
multiplied  by  the  number  of  calibrations  to  be  completed,  it  was  apparent 
that  too  much  time  would  be  spent  doing  calibrations. 

All  that  is  needed  from  the  calibrations  are  frequency  responses  from  various 
components,  the  only  exception  being  crystal  diode  detectors  where  aDplied 
power  variation  is  also  needed.  Obtaining  these  variations  is  simple  for  a 
program  that  loops  on  the  frequency  and/or  power  and  measures  the  response  of 
the  component  under  test.  This  is  the  basis  for  all  the  calibration  acquisi¬ 
tion  software. 

The  problem  with  the  above  acquisition  scheme  is  the  frequency  response  of 
the  various  monitor  components  (i.e.,  power  sensors  and  directional  couplers) 
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to  microwaves.  The  solution  to  this  problem  is  to  have  the  reduction  soft¬ 
ware  read  the  raw  data  files  and  remove  the  response  of  the  power  sensors  and 
any  directional  couplers  used  to  obtain  the  raw  data.  This  process  gives  a 
better  method  of  knowing  what  signal  had  to  be  applied  to  a  component  in 
order  to  get  the  measured  data. 

A  side  benefit  from  completing  our  calibrations  with  this  hardware  is  por¬ 
tability.  The  calibration  cart  can  be  used  anywhere  in  the  laboratory  area 
and  calibrate  devices  in  place.  This  provides  better  calibrations  and  also 
allows  for  modifications  of  experiments  without  having  to  completely 
disassemble  and  rebuild  the  experiment  with  independently  calibrated 
hardware. 

HARDWARE 

The  hardware  used  to  obtain  the  calibration  data  works  under  the  premise 
that,  if  enough  data  points  are  used,  one  can  adequately  approximate  the  real 
response  function.  The  hardware  (Fig.  1)  necessary  to  accomplish  this  tabu¬ 
lation  of  data  includes:  (1)  a  controller  (Z100  microcomputer)  capable  of 
large  amounts  of  data  storage;  (2)  an  oscillator  (HP  835QA)  with  variable 
power  and  frequency;  (3)  various  power  sensors  (HP  8478B,  848 1A,  8484A)  to 
measure  the  applied  power  to  the  device  under  test;  (4)  an  interface  [a  BIRA 
analog-to-digital  converter  (ADC)  housed  in  a  CAMAC  crate]  to  transfer  the 
data  from  the  power  sensors  to  the  storage  device;  and  (5)  a  computer  (VAX 
11/730)  to  take  all  the  various  sized  files  and  unfold  the  imbedded  sensor 
response  funtions. 

The  basic  change  to  the  above  list  is  the  controller.  Previously  a  Tektronix 
4051  microcomputer,  with  GPIB  and  RS-232C  interfaces,  was  used.  We  are  now 
using  a  Zenith  Z-100  microcomputer  with  a  built-in  RS-232C  interface  and  an 
add-on  National  Instruments  GPIB  interface.  The  software  driver  was 
purchased  from  National  Instruments  for  ZBASIC  running  under  the  MS-DOS 
operating  system. 

The  reason  for  using  the  Z100  instead  of  the  4051  was  its  increased  storage 
ability.  The  Z100  is  equipped  with  a  10-megabyte  Winchester  drive  and  a 
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Figure  1.  Z100  calibration  system 


360-kilobyte  floppy  drive.  The  programs  are  maintained  on  the  Winchester 
disk,  and  the  date  files  are  kept  on  various  floppy  diskettes  for  crystal, 
attenuators,  cables,  and  other  designated  diagnostics. 

The  one  shortfall  of  the  Z100  is  the  limitation  in  memory  space  for  data 
storage.  The  amount  of  data  required  by  the  reduction  software  typically 
exceeds  the  64-kilobyte  limit  and  as  a  result,  our  reduction  software  is  run 
on  the  VAX  11/750.  Having  to  use  the  VAX  for  data  reduction  requires  extra 
time  when  the  calibration  cart  must  move  about  the  lab.  However,  since  many 
IBM-AT  type  computers  are  available  (with  terminal  software),  one  can  get  to 
the  VAX  by  taking  the  floppy  disk  to  any  of  the  microcomputers. 

SOFTWARE 

Z100 

The  software  for  the  Z100  is  essentially  the  same  as  that  used  on  the  4051, 
modulo  language  differences.  We  have  moved  the  software  for  the  crystal 
calibrations,  directional  coupler  calibrations,  and  "attenuation"  calibra¬ 
tions.  Hookup  diagrams  are  shown  in  Figs.  2,  3,  and  4. 

Figure  2  indicates  the  intended  setup  for  the  program  DCCAL.BAS.  The 
through-power  sensor’s  recorder  outDut  is  hooked  to  port  1  of  the  BIRA  ADC 
bus.  The  pickoff-power  sensor’s  recorder  output  goes  to  port  2  of  the  BIRA 
ADC  bus. 

Figure  3  indicates  the  hookups  for  measurements  made  using  the  program 
DB.BAS.  The  incident -power  sensor  is  attached  to  port  1  of  the  BIRA  ADC,  the 
reflected-power  sensor  is  attached  to  port  2  of  the  BIRA  ADC,  and  the 
transmitted-power  sensor  is  attached  to  port  3  of  the  BIRA  ADC. 

Figure  4  indicates  the  proper  set  up  for  using  XTAL.BAS.  The  incident  and 
reflected  sensors  are  connected  as  for  DB.BAS.  The  crystal  must  be 
calibrated  by  driving  the  same  load  it  sees  during  actual  experiments.  In 
this  work,  that  is  50  12.  An  artifact  of  driving  a  50-£2  load  is  to  cause  the 
voltage  level  to  be  between  0  and  -100  mV.  This  voltage  is  to  be  measured  by 
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an  ADC  with  12  bits  of  resolution  for  a  0-  to  +5-V  range.  To  get  around 
using  only  the  low  bits  of  the  ADC  and  the  polarity  problem,  a  simple 
inverting  amplifier  circuit  (using  a  741  operational  amplifier)  is  used.  The 
gain  is  set  to  ~  58  (actually  35.3  dB  ±  0.2  dB  =  20.  *  log  IV  tAinU  over 
the  output  limits  of  the  sweep  oscillator.  The  gain  is  then  divided  out  in 
the  program  and  the  quotient  is  stored  as  the  measured  value. 

These  proarams  have  considerable  flexibility  built  into  them.  An  example  is 
shown  in  Fig.  5.  This  setup  was  used  to  calibrate  cables  outside  the  fre¬ 
quency  range  of  the  directional  couplers.  Balanced  power  splitters  are 
available  which  may  be  used  to  DC  frequencies.  The  splitter  was  first 
calibrated  and  found  to  have  0  dB  (±  0.2  dB)  relative  attenuation  between 
outputs;  this  meant  the  splitter  could  obtain  directly  the  frequency  response 
of  cables,  or  any  other  component  response  down  to  10  MHz  (the  limit  of  the 
sweep  oscillator). 

Program  listings  are  included  in  the  appendix.  As  noted  in  Ref.  1,  segments 
in  the  programs  can  change  the  parameter  range  of  the  calibration.  The  pararm 
eters  one  may  want  to  change  ares  (1)  power  level  or  range  of  power  levels 
applied,  (2)  frequency  range  covered,  (3)  number  of  data  points  measured  both 
in  freauency  and  power,  and  (4)  power  meter  scale  factors.  In  the  Drogram 
DCCAl.BAS,  line  1130  specifies  the  start  frequency  (FI)  and  stop  frequency 
(F2) ,  the  applied  power  level  (PO),  and  the  number  of  data  points  (NO). 

Also,  in  line  1227  the  power  meter  scale  factors  are  preset.  Parameter  M(l) 
is  the  scale  factor  needed  to  convert  the  voltage  output  read  on  ADC  port  1 
to  the  indicated  power  level  on  the  meter,  and  M(2)  is  the  scale  factor  for 
ADC  port  2. 

The  scale  factors  are  determined  in  the  following  manner.  The  power  meter 
manuals  state  that  the  recorder  output  is  a  DC  voltage  between  0  and  1.2  V 
regardless  of  power  meter  scale.  For  example,  a  voltage  of  0.9  V  is  output 
for  90  yW  on  the  100-yW  scale,  900  yW  on  the  1-mW  scale,  and  9  mW  on  the 
10-mW  scale.  The  scale  factor  is  that  number  which,  when  multiplied  with  the 
voltage  from  the  recorder  output,  gives  a  power  level  in  mW.  In  the  case  of 
90  yW,  the  scale  factor  would  be  0.1;  for  900  yW  the  scale  factor  would  be  1* 
and  for  9  mW  the  scale  factor  would  be  10. 
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The  parameter  block  in  ATTN. BAS  is  found  at  lines  1375  to  1430.  The  param¬ 
eter  CO  is  a  counter  used  for  delay  loops  and  can  be  ignored.  The  parameter 
block  for  XTAL.BAS  is  found  at  lines  1490  to  1580.  The  additional  parameters 
are  PI  and  P2,  the  start  and  stop  power  level,  N1  the  number  of  power  levels 
to  take  data,  YMAX  a  typical  maximum  voltage  output  by  the  crystal  diodes, 
and  G1  the  gain  of  the  amplifier  box.  All  other  parameters  are  requested  by 
the  program  and  need  no  explanation. 

VAX 


The  programs  for  the  VAX  are  also  included  in  the  appendix.  Notice  that  for 
the  files  accessed  by  programs  OB. FOR  and  XTAL.FOR,  one  can  have  applied 
power  computed  in  two  ways:  (1)  using  only  the  measured  incident  power;  and 
(2)  using  the  difference  of  the  incident  and  reflected  power.  Case  [1)  is 
more  useful  for  pulsed  microwave  measurements.  Case  [2)  is  more  useful  in 
CW  microwave  measurements. 

These  programs  use  the  subroutines  DCCOR  and  TCOR  to  remove  the  frequency 
response  of  the  power  sensors  and  calibrated  directional  couplers  used  in  the 
data  acquisition  to  calculate  the  applied,  reflected,  or  transmitted  power 
signal.  These  subroutines  work  from  data  files  and  must  interpolate  between 
the  data  points.  The  data  is  interpolated  by  a  cubic  spline  method  to  give 
better  accuracy  than  a  simple  linear  interpolation. 

The  program  DRC.FOR  uses  the  output  file  from  DCCAL.BAS  to  give  frequency 
response  of  directional  couplers,  or  cables.  The  output  of  DRC.FOR  is  a  file 
DRATTN.Pat  which  should  be  renamed  to  indicate  the  calibrated  device. 


The  program  DB.FOR  uses  data  files  created  by  ATTN. BAS.  The  reduced  data  are 
saved  in  a  file  called  SAMPLE.DAT.  This  file  should  also  be  renamed  for 
later  use. 

The  programs  XTAL.FOR  and  PSWP.FOR  are  used  to  get  calibration  data  from 
XTAL.BAS.  XTAL.FOR  reads  the  raw  data  file  from  XTAL.BAS  and  folds  out  the 
frequency  response  of  the  directional  couplers  and  power  sensors  to  get  back 
to  the  actual  power  applied  to  the  crystal.  After  one  decides  how  the 
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applied  power  is  to  be  calculated,  the  files  XTALV.DAT  and  XTALP.DAT  are 
created.  These  files  contain  information  on  the  voltage  and  power  measured 
in  matrix  form.  The  columns  of  each  file  are  “constant"  power  applied,  and 
the  rows  are  “constant"  frequency  applied.  These  files  should  be  renames  for 
later  use.  The  renamed  files  are  used  by  PSWP.FOR  to  give  a  power  sweep  at  a 
chosen  frequency.  If  the  chosen  frequency  is  not  one  of  the  measured  fre¬ 
quencies,  the  data  in  both  files  are  interpolated  by  computing  a  cubic  spline 
on  the  frequencies  immediately  above  and  below  and  a  computed  set  of  numbers 
are  generated  to  fill  in  XTALV  and  XT ALP.  These  computed  numbers  are  then 
used  to  give  back  to  the  user  the  power  level  necessary  to  generate  the 
measured  voltage. 

In  running  all  the  data  reduction  software,  graphs  are  generated  at  many 
places  that  are  useful  for  later  quick  calculations  or  when  looking  for  “bad" 
data.  A  hardcopy  might  be  helpful.  The  programs  also  create  data  files 
which  are  tables  of  the  data  used  to  draw  the  graphs.  These  files  are  useful 
when  more  accurate  information  is  needed  for  calculations. 

CONCLUSIONS 

This  report  has  discussed  an  improved  microwave  calibration  system.  The 
three  areas  covered  were*  (1)  the  hardware  used  to  build  the  calibration 
system;  (2)  the  software  used  to  acquire  the  raw  calibration  data;  and  (3) 
the  software  used  to  reduce  the  raw  data  and  obtain  the  "true**  calibration 
response  functions  of  our  microwave  detectors. 

The  system  does  operate  as  advertised  on  the  date  of  this  report.  Modifica¬ 
tions  to  the  system  may  occur  in  the  future,  and  will  be  properly  documented 
at  that  time. 
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DCCAL.BAS 


APPENDIX 

PROGRAM  LISTINGS 


10  CLEAR  ,60000!  '  BASIC  Declarations 

20  IBINIT1  =  60000! 

30  IBINIT2  =  IBINIT1  +  3 

40  BLOAD  *  bib . m* , IBINIT 1 

50  CALL  IBINIT1 (IBFIND , IBTRG, IBCLR, IBPCT, IBSIC  , 

IBLOC , IBPPC , IBBNA , IBONL , IBRSC . 

IBSRE . IBRSV, IBPAD , IBS AD , I  BIST , 

IBDMA, I BEOS , IBTMO , IBEOT , IBRDF , IBWRTF) 
60  CALL  IBINIT2 (IBGTS , IBCAC , IB WAIT , IBPOKE , IBWRT , 

IBWRTA, IBCMD , IBCMDA, IBRD , IBRDA, 

IBSTOP , IBRPP , IBRSP , IBDI AG , IBXTRC , 

IBRD  I  ,  IBWRT  I  ,  IBRDIA,  IBWRTIA,  IBSTA7. , 
IBERR%. IBCNTX) 

1000  CLS 

1010  WRT*  = ‘ SWEEPER ' : CALL  IBFIND (WRT* , SWEEPERX) 

1020  WRT*  =  ‘  CAMAC  '  :  CALL  IBFIND  (WRT*  ,  CAMAC7.) 

1030  WRT*  = ‘ GPIBO ' : CALL  IBFIND ( WRT* , GPIB0%) 

1040  CALL  IBCLR  (SWEEPER*/.) 

1045  REM  TURN  RF  OFF  AND 

1046  REM  TURN  CW  FILTER  AND  INTERNAL  LEVELING  ON 
1050  WRT** * RFO *: CALL  IBWRT  (SWEEPER?. ,  WRT*) 

1052  WRT*=  * A1 * : CALL  IBWRT (SWEEPER% , WRT*) 

1054  WRT**  *  FL1  *:  CALL  IBWRT (SWEEPERS , WRT*) 

1060  REM  HEADER  BLOCK  ENTRY 
1070  PRINT  ‘ENTER  THE  FILENAME* 

1080  INPUT  V* 

1090  OPEN  *0* ,*1 ,V* 

1100  PRINT  ‘ENTER  TEST  LABEL  (IN  SINGLE  QUOTES) * : INPUT  V* 
1110  10=0 : WRITE  *1 ,10 
1120  DIM  D (3) ,M(2) 

1125  REM  DEFINE  FREQUENCY  RANGE,  POWER  LEVEL, 

1126  REM  AND  NUMBER  OF  DATA  POINTS 
1130  F1=.5:F2=12! :P0=11! :N0=30 

1140  PRINT  *1  ,V*;  ,- 'POWER  LEVEL  (DBM)’";P0 
1145  REM  SET  FREQUENCY  STEP  SIZE 
1150  LET  SO* (F2-F1) / (NO- 1) 

1160  PRINT  ’ENTER  ID  FOR  DIRECTIONAL  COUPLER' : INPUT  D(3) 
1170  PRINT  'ENTER  ID  FOR  THROUGH  THERM  MOUNT': INPUT  T1 
1180  PRINT  'ENTER  ID  FOR  PICKOFF  THERM  MOUNT': INPUT  T2 
1190  REM  SEND  HEADER  TO  11/60 

1200  PRINT  *1 ,' ’DIRECTIONAL  COUPLER  ID*’*:WRITE  *1,D(3) 

1210  PRINT  *1,* 'THROUGH  THERMISTOR  MOUNT  ID»’*:WRITE  *1,T1 
1220  PRINT  *  1  , " ’ PICKOFF  THERMISTOR  MOUNT  ID*’ *: WRITE  «1,T2 
1225  REM  PRESET  POWER  METER  SCALE  FACTORS 
1227  M(l) =10:M(2) =1  ! 

1230  PRINT  'THE  FOLLOWING  SETUP  IS  USED' 

1240  PRINT  'THE  THROUGH  POWER  IS  MONITORED  ON  PORT  1" 

1250  PRINT  'THE  PICKOFF  POWER  METER  IS  ON  A  SCALE  OF  *;M(1) 
1260  PRINT  'THE  PICKOFF  POWER  METER  IS  ON  A  SCALE  OF  *;M(2) 
1270  PRINT  'ENTER  A  'Y*  WHEN  SET' 
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1280 

1290 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1485 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1580 

1590 

1690 

1700 

1710 

1711 

1712 

1713 

1714 
1720 
1780 
1785 
1790 
1800 
1810 
1820 
1830 
1835 
1840 
1850 
1860 
1865 
1870 
1880 
1881 
1890 


INPUT  V* 

IF  V*=CHR*(89)  OR  V*=CHR* ( 1 2 1 )  THEN  1310  ELSE  1280 
PRINT  *1 , ‘ ’FREQ(HZ)  THROUGH  POWER(W)  PICKCFF  POWER(W)’ 
CALL  I BC  LR(S  WEEPER /() 

WRT#=*FA*  :CALL  IBWRT  (SWEEPER?. ,  WRT*) 

WRT*=STR*(F1) : CALL  IBWRT (SWEEPER% , WRT*) 

WRT*  = *  GZ  *  : CALL  IBWRT (SWEEPER* , WRT*) 

WRT*  = ' FB " : CALL  I BWRT (SWEEPER% , WRT* ) 

WRT*=STR*  (F2)  :  CALL  IBWRT  ( SWEEPER?. ,  WRT*) 

WRT*  =  "GZ" : CALL  IBWRT (SWEEPER% , WRT* ) 

WRT*  =  *  PL ‘ : CALL  I BWRT ( SWEEPER* , WRT* ) 

WRT*=STRS  (PO)  :  CALL  IBWRT  ( SWEEPER*/. ,  WRT* ) 

WRT*  =  "DM"  :  CALL  IBWRT  (SWEEPER*/. ,  WRT*) 

WRT*=  '  SM"  :  CALL  IBWRT  (SWEEPER*/. ,  WRT*) 

WRT*  =  STR*  (FI)  :  CALL  I  BWRT  ( SWEEPER*/. ,  WRT* ) 

WRT*=  *  GZ ‘ : CALL  I BWRT ( SWEEPER* , WRT* ) 

WRT*=  *  SS ' : CALL  IBWRT (SWEEPER% , WRT*) 

WRT*  =  STR*  ( 1000  !  *S0)  :  CALL  IBWRT  (SWEEPER*/. ,  WRT*) 

WRT*= " MZ ’ : CALL  IBWRT (SWEEPER* , WRT*) 

CLS 
N2=  1 

REM  SET  UP  GRAPH  AXES 
LINE  (25 , 10) - (25 ,210) 

LINE  -(565,210) 

CALL  IBCLR (SWEEPER*) 

WRT*  =  '  RF 1  *  :  CALL  IBWRT  (SWEEPER?. ,  WRT*) 

FOR  1=1  TO  NO 
F=F 1 + ( I - 1 ) *S0 
FOR  J=1  TO  2 
A=J 

Q=M(J) *5/4095 
FOR  V=1  TO  N2 
NEXT  V 
GOSUB  1890 

IF  Y/M( J) >1.2  THEN  1713 
IF  Y/M( J) > . 099  THEN  1720 
GOSUB  2200 
GOTO  1580 
D  ( J) =Y 

YMAX= .9* ( 10* (P0/10) -1) 

REM  PLOT  DATA  POINTS 

PSET  ( ( (F-Fl) / (F2-F1) ) *540+25 , 2 10* ( 1-D (J) / YMAX) ) 

NEXT  J 

WRITE  *l,F«10“9,D(l)»10"-3,D(2)»10"-3 
CALL  IBCLR  (SWEEPER*/.) 

WRT* »* SM* : CALL  IBWRT (SWEEPER* , WRT*) 

WRT*=  *  UP  * : CALL  IBWRT (SWEEPER% , WRT*) 

NEXT  I 

REM  CALL  IBCLR (SWEEPER%) 

WRT*=  *  RFO  '  :  CALL  IBWRT (SWEEPER* , WRT*) 

CALL  IBLOC (SWEEPER*) 

WRITE  *1,-999,-999,-999,-999 

CLOSE  *1 

END 

REM  CAMAC  DATA  FETCH 
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1900  REM  CALL  IBSIC (GPIBO*) 

1910  REM  WRT*  = ' C AMAC ‘ : REM  CALL  IBFIND ( WRT* , CAMAC*) 

1920  WRT*=CHR* (30) +  CHR* (0) +  CHR* ( 17) +  CHR* (0) +CHR* (0) +CHR* (0) 
1930  CALL  IBWRT (CAMAC* , WRT*) 

1940  WRT*=CHR* (14) +  CHR* (0) +  CHR* (26) +  CHR* (0) +CHRS (0 ) +CHR* ( A) 
1950  CALL  I BWRT ( CAMAC* , WRT* ) 

1960  WRT*=CHR* (14) +CHR* (0) +  CHR* (16) +CHR* (0) +  CHR* (0) +  CHR* (A) 
1970  CALL  IBWRT (CAMAC* ,  WRT*) 

1980  WRT*=CHR* (14) +CHR* (0) +CHR* (0) 

1990  CALL  IBWRT  (CAMAC* , WRT*) 

2000  CNT*=3 : CALL  IBRDI (CAMAC* ,C* (0) .CNT*) 

2010  B1*=C*(0) /256 

2020  Y=Q* (B 1**256+ ABS (C% ( 1) ) +3) 

2030  WRT*=CHR* ( 14) +  CHRS (0) +  CHR* (24) +  CHR* (0) +  CHRS (0) +  CHRS (A) 
2040  CALL  IBWRT  (CAMAC* , WRT* ) 

2050  WRT*=CHR* (30) +CHR* (0) +  CHR* ( 17) +  CHR# ( 0 ) +  CHR* (0) +  CHR* (96) 
2060  CALL  I BWRT (CAMAC* , WRT* ) 

2070  IF  Y<0  THEN  1890 
2080  RETURN 
2090  REM  METER  ZERO 
2100  CALL  IBCLR (SWEEPER*) 

2110  WRT*  = ‘ *  FO ’ : CALL  IBWRT (SWEEPER* , WRT*) 

2120  LOCATE  1 , 1 

2130  PRINT  ’ENTER  A  Y  WHEN  METERS  ZEROED INPUT  V* 

2140  IF  V* <  > ' Y *  OR  V*<  > "y '  THEN  2190 
2150  WRT*  = ’ RF 1 ’ : CALL  IBWRT (SWEEPER* , WRT*) 

2160  FOR  V=1  TO  25 
2170  NEXT  V 
2200  REM  SCALE  CHANGE 
2210  LOCATE  1 , 1 

2220  PRINT  ’SCALE  SETTING  FOR  PORT  ’;J;*  =  ’;M(J) 

2230  PRINT  ’ENTER  NEW  SCALE  VALUE’: INPUT  M(J) 

2240  RETURN 
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1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1085 

1086 
1090 
1092 
1094 
1  100 
1  1  10 
1120 
1  130 
1140 
1150 
1160 
1170 
1 180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 


CLEAR  ,60000!  ’  BASIC  Declarations 

IBINIT1  =  60000! 

IBINIT2  =  IBINIT1  +  3 
BLOAD  ’ bib . m’ .IBINITl 

CALL  IBINITl (IBFIND , IBTRG , IBCLR , IBPCT, IBSIC , IBLOC , 
IBPPC , IBBNA, IBONL, IBRSC, IBSRE, IBRSV, 
IBPAD , IBS AD , IBIST , IBDMA, IBEOS , IBTMO , 
IBEOT , IBRDF , IBWRTF) 

CALL  IBINIT2 ( IBGTS , IBCAC , IBWAIT , IBPOKE , IBWRT , IBWRTA , 
IBCMD , IBCMDA, IBRD , I BRDA , I BSTOP , IBRPP , 
IBRSP , IBDIAG , IBXTRC , IBRD I , IBWRT I , IBRDIA, 
IBWRTIA,  I BSTA%  ,  IBERR7. ,  IBCNT7.) 

CLS 

REM  THE  DETECTOR  HOOK-UP  IS  AS  FOLLOWS 
REM  PORT  1  IS  THE  INCIDENT  POWER  MONITOR 
REM  PORT  2  IS  THE  REFLECTED  POWER  MONITOR 
REM  PORT  3  IS  THE  TRANSMITTED  POWER  MONITOR 
WRT*  =  ‘  SWEEPER’  :  CALL  IBFIND  ( WRT*  ,  SWEEPER"/.) 

WRT  $  =  ‘  C  A  MAC  *  :  CALL  IBFIND  (WRT*  .CAMAC*/.) 

WRT*  =  “ GPIBO  *  :  CALL  IBFIND  (WRT*  .GPIBOTi) 

CALL  I BCLR ( SWEEPERS ) 

REM  TURN  RF  SOURCE  OFF  AND  TURN  ON  CW  FILTER 

REM  AND  INTERNAL  LEVELING 

WRT*  =  * RFO ' : CALL  IBWRT (SWEEPER* , WRT*) 

WRT*  =  * A1 * : CALL  IBWRT (SWEEPERX , WRT*) 

WRT*  = ' FL 1 " : CALL  IBWRT (SWEEPER* , WRT*) 

REM  HEADER  BLOCK  ENTRY 
PRINT  *  ENTER  THE  FILENAME* 

INPUT  V* 

OPEN  *0* ,*1 . V* 

PRINT  *  ENTER  TEST  LABEL  (IN  SINGLE  QUOTES)  *: INPUT  V* 

10  =  0 

PRINT  *  ENTER  SAMPLE  ID*  * : INPUT  LI 
PRINT  ‘ENTER  ID  FOR  INC  DIR  COUPLER *: INPUT  D1 

PRINT  "ENTER  ID  FOR  REF  DIR  COUPLER* : INPUT  D2 

PRINT  *  ENTER  ID  FOR  INC  THERM  MOUNT’: INPUT  T1 

PRINT  *  ENTER  ID  FOR  REF  THERM  MOUNT*: INPUT  T2 

PRINT  *  ENTER  ID  FOR  TRANS  THERM  MOUNT’: INPUT  T3 
REM  SEND  HEADER  TO  11/60 
PRINT  *i ,V* 

WRITE  *1,10 

PRINT  *1 , ’ ’CABLE  ID’ ’ 

WRITE  *  1 , L 1 

PRINT  *  1  ,  ’  ’  INC  DIR  COUPLER  ID” 

WRITE  *1 ,D1 

PRINT  *  1  ,  ”  REF  DIR  COUPLER  ID” 

WRITE  *1 ,D2 

PRINT  *  1  ,  ’  ’  INC  THERM  ID” 

WRITE  *  1 ,T1 

PRINT  *  1  ,  ”  REF  THERM  ID” 

'WRITE  *  1 ,T2 

PRINT  *1,’  ’TRANS  THERM  MOUNT  ID” 
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1360 

1370 

1375 

1380 

1390 

1400 

1405 

1410 

1430 

1450 

1455 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1562 

1563 

1564 

1565 

1566 

1567 
1570 
1580 
1590 
1600 
1610 
1615 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1700 
1710 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1940 
1945 
1950 
1960 


WRITE  #1 ,T3 

PRINT  *  1 , * ’ FREQ (HZ)  INC  POW(W)  REF  POW(W)  TRANS  POW(W) 
REM  SET  CALIBRATION  PARAMETERS 
FI*  1  i 
F2= 18 ■ 

P0=  13 

REM  NUMBER  OF  FREQUENCY  STEPS 

N0  =  200 

C0=1 

DIM  D (3) , M( 3 ) 

REM  SET  FREQUENCY  STEP  SIZE 
S0= (F2-F1) / (NO-1) 

CALL  IBCLR  (SWEEPER*/.) 

WRT*=*FA*  :  CALL  IBWRT (SWEEPER* , WRT*) 

WRT*=STR*  (FI):  CALL  I  BWRT  ( SWEEPER*/. ,  WRT* ) 

WRT*  =  *  GZ  *  :  CALL  I  BWRT  (SWEEPER*/. ,  WRT*) 

WRT*=  * FB " : CALL  I BWRT (SWEEPER* , WRT* ) 

WRT*  =  STR*  ( F2 )  :  CALL  I  BWRT  (SWEEPER*/. .  WRT*) 

WRT*=  "  GZ  ’  :  CALL  I  BWRT  (SWEEPER*/. ,  WRT* ) 

WRT*=  *  PL  *  :  CALL  I  BWRT  (SWEEPER*/. ,  WRT*) 

WRT*  =  STR*  (PO)  :  CALL  IBWRT  ( SWEEPER?. ,  WRT* ) 

WRT*=*DM* :CALL  IBWRT (SWEEPER% , WRT*) 

WRT*=  *  SM* :  CALL  IBWRT (SWEEPER* . WRT*) 

WRT*  =  STR*  (FI)  :  CALL  IBWRT  ( SWEEPER*/. ,  WRT* ) 

WRT*  =  *  GZ  *  :  CALL  IBWRT (SWEEPER* , WRT*) 

WRT*=  * SS * :  CALL  IBWRT (SWEEPER* , WRT*) 

WRT*=STR* (S0« 1000) :  CALL  IBWRT (SWEEPER* , WRT*) 

WRT* *  *MZ* :  CALL  IBWRT (SWEEPER* .WRT*) 

PRINT  "ENTER  THE  SCALE  FACTORS  FOR  THE  POWER  METERS’ 


PRINT  ’ SCALE  FACTOR  FOR  PORT  1 *: INPUT  M(l) 

PRINT  ’SCALE  FACTOR  FOR  PORT  2*: INPUT  M(2) 

PRINT  ‘SCALE  FACTOR  FOR  PORT  3*: INPUT  M(3) 

CLS 

REM  DRAW  AXES 

LINE  (25 , 10) - (25 ,210) 

LINE  -(565,210) 

CALL  IBCLR (SWEEPER*) 

WRT*  =  * RF 1  * : CALL  IBWRT (SWEEPER* , WRT*) 

FOR  1=1  TO  NO 
F  =  F  1  + { I -  1 ) *S0 
FOR  K= 1  TO  3 
A  =  K 

Q=M(K) *5/4095 
FOR  V=1  TO  CO 
NEXT  V 
GOSUB  2060 

IF  Y/M(K) >1.2  THEN  1860 
IF  Y/M(K) > . 099  THEN  1880 
GOSUB  2260 
GOTO  1700 
D (K) = Y 

YMAX= (10* (PO/ 10) ) -  1 

REM  PLOT  DATA  POINT  ON  GRAPH 

PSET  ( ( (F-Fl) / (F2-F1) ) *540+25 ,210* ( l-D(K) / YMAX) ) 
NEXT  K 
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1970  WRITE  *1  ,F«  10*9  ,D ( 1) *  10*-3 ,D(2) *  10*-3 ,D (3) « 10*-3 
1980  CALL  IBCLR(SWEEPER%) 

1990  WRT*  = ‘ SM" : CALL  I BWRT ( SWEEPER* . WRT* ) 

2000  WRT*= 'UP' :CALL  I  BWRT  (SWEEPER*/. ,  WRT*) 

2010  NEXT  I 

2030  WRT*='RF0' :CALL  IBWRT  ( SWEEPER'/. ,  WRT* ) 

2035  CALL  IBLOC (SWEEPER%) 

2040  WRITE  #1,-999,-999,-999,-999 

2050  CLOSE  #1 

2051  END 

2060  REM  CAMAC  DATA  FETCH 
2070  REM  CALL  IBSIC  (GPIB07.) 

2080  REM  WRT*= ' CAMAC ': REM  CALL  IBFIND (WRT* , CAMAC%) 

2090  WRT*=CHR* (30) +  CHR* (0) +  CHR* (17) +  CHR* (0) +  CHR* (0) +CHR* (0) 
2100  CALL  IBWRT  (CAMAC'/. ,  WRT*) 

2110  WRT*=CHR* (14) +  CHR* (0) +CHR* (26) +  CHR* ( 0 ) +CHR* (0) +  CHR* (A) 
2120  CALL  IBWRT (CAMACX , WRT*) 

2130  WRT*=CHR* (14) +CHR* (0) +  CHR* (16) +CHR* (0) +CHR* (0) +CHR* (A) 
2140  CALL  I  BWRT  (CAMAC*/.,  WRT*) 

2150  WRT*=CHR* (14) +CHR* (0) +CHR* (0) 

2160  CALL  IBWRT  (CAMAC3C ,  WRT* ) 

2170  CNT%=3 : CALL  IBRD I (CAMACX , C% ( 0 ) , CNT%) 

2180  B1/C  =  C%  (0)  /256 

2190  Y=Q* (B1X«256+ABS (C% ( 1 ) ) +3) 

2200  WRT*=CHR* ( 14) +CHR* (0) +  CHR* (24) +  CHR* (0) +CHR* (0) +  CHR* (A) 
2210  CALL  IBWRT (CAMACX , WRT*) 

2220  WRT*=CHR* (30) +CHR* (0) +  CHR* ( 17) +  CHR* (0) +CHR* (0) +CHR* (96) 
2230  CALL  IBWRT  (CAMAC7.  ,  WRT*) 

2240  IF  Y<0  THEN  2060 
2250  RETURN 
2260  REM  SCALE  CHANGE 
2270  LOCATE  1 , 1 

2280  PRINT  'SCALE  SETTING  FOR  PORT*;K;'  =  *;M(K) 

2290  PRINT  'ENTER  NEW  SCALE  VALUE': INPUT  MOO 

2300  RETURN 

2310  REM  METER  ZERO 

2320  CALL  IBCLR (SWEEPERX) 

2330  WRT*='RF0" :CALL  I BWRT ( SWEEPERX , WRT* ) 

2340  LOCATE  1,1 

2350  PRINT  'ENTER  A  Y  WHEN  METERS  ZEROED ': I NPUT  V* 

2360  IF  V*<  > ' Y'  OR  V*<>'y'  THEN  2410 
2370  WRT*=  '  RF1 '  :  CALL  I  BWRT  ( SWEEPER'/. ,  WRT* ) 

2380  FOR  V=1  TO  25 
2390  NEXT  V 
2460  RETURN 
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XTAL. BAS 


1000 

1010 

1020 

1030 

1040 

1050 

1060 

1065 

1066 
1070 
1080 
1090 
1  100 
1  1 10 
1120 
1  130 
1140 
1150 
1  160 
1  170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1380 
1390 


CLEAR  ,60000!  '  BASIC  Declarations 

IBINIT1  =  60000! 

IBINIT2  =  IBINIT1  +  3 
BLOAD  ' bib . m* , IBINIT1 

CALL  IBINIT1 (IBFIND , IBTRG, IBCLR , IBPCT , IBSIC , I BLOC , 
IBPPC , IBBNA, IBONL, IBRSC , IBSRE , IBRSV, 
IBPAD , IBSAD , I  BIST , IBDMA , IBEOS , IBTMO , 
IBEOT , IBRDF , IBWRTF) 

CALL  IBINIT2 ( IBGTS , IBCAC , IBWAIT , IBPOKE , IBWRT , IBWRTA , 
IBCMD , IBCMDA , IBRD , IBRDA , IBSTOP . IBRPP , 
IBRSP , IBDIAG. IBXTRC . IBRD I , IBWRT I , IBRD  I A , 
IBWRTIA,  IBSTA7. .  IBERR7. ,  IBCNT*/.) 

THIS  PROGRAM  RUNS  A  DOUBLE  FOR-NEXT  LOOP 
THE  OUTER  LOOP  IS  A  POWER  SWEEP 
THE  INNER  LOOP  IS  A  FREQ  SWEEP 


IBPAD , 
IBEOT 
IBINIT2 (IBGTS 
IBCMD 
IBRSP 


REM  THIS  PROGRAM  R 
REM  THE  OUTER  LOOP 
REM  THE  INNER  LOOP 


BD*= ’ SWEEPER" : CALL  IBFIND (BD* , SWEEPER%) 

BD*  = ‘ GPIBO ' : CALL  IBFIND (BD* .GPIBO) 

CALL  IBCLR  (SWEEPER*/.) 

REM  TURN  RF  OFF  AND  TURN  ON  INTERNAL  LEVELING 
REM  AND  CW  FILTER 

WRT*  = ’ RFO ‘ : CALL  IBWRT ( SWEEPERS , WRT* ) 

WRT*  =  *  A1  ’  :  CALL  IBWRT  (SWEEPER?. .  WRT*) 

WRT*  = ‘ FL 1  * : CALL  IBWRT (SWEEPER* , WRT*) 

PRINT  ’ENTER  THE  FILE  NAME- 
INPUT  V* 

OPEN  -0- ,*1,VI 

PRINT  - ENTER  THE  TEST  LABEL ( IN  SINGLE  QUOTES) ’ 

INPUT  V# 

PRINT  *1  ,V* 

REM  SEND  DATA  FLAG 

WRITE  *1 ,0 

REM  HEADER  BLOCK 

DIM  D  (3)  ,  M  ( 3)  ,  T  ( 2 )  ,C27.(3) 

PRINT  'ENTER  XTAL  ID*-: INPUT  D { 1 ) 

PRINT  'ENTER  THE  INC  DIRECTIONAL  COUPLER  ID*-: INPUT  D(2) 

PRINT  - ENTER  THE  REF  DIRECTIONAL  COUPLER  ID*-: INPUT  D(3) 

PRINT  ’ENTER  THE  INC  THERM  MOUNT  ID*-: INPUT  T ( 1 ) 

PRINT  ‘ ENTER  THE  REF  THERM  MOUNT  ID*-: INPUT  T(2) 

PRINT  *1 , - ’XTAL  ID* ’ - 
WRITE  #1  ,D(  1) 

PRINT  *1 , - ’ INC  DRC  ID* ’ * 

WRITE  *1  , D ( 2) 

PRINT  *1 , * ’REF  DRC  ID* ’ - 
WRITE  *1  ,  D  ( 3 ) 

PRINT  *  1 , -  ’ INC  THERM  ID* ’ - 


: INPUT  T (2) 


PRINT  *  1 , -  ’ INC  THERM  ID* ’ - 
WRITE  *1 , T ( 1 ) 

PRINT  *  1 , - ’ REF  THERM  ID*’- 
WRITE  *1  , T ( 2 ) 

REM  POWER  METER  SCALE  FACTORS 
M(  3 )  =  1 

PRINT  ‘ENTER  A  Y  IF  THE  XTAL  GIVES 


POSITIVE 
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1391 

1395 

1396 

1397 
1400 
1420 
1430 
1440 
1450 

1460 

1461 
1470 
1480 
1490 
1495 
1500 
1510 
1515 
1520 
1530 
1535 
1540 
1545 
1550 
1560 
1570 
1575 
1580 
1585 
1590 
1600 
1605 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
T790 
1800 
1805 
1810 
1820 


PRINT  ’SIGNAL  INTO  50  OHMS’ 

PRINT  ’MAKE  SURE  THE  GAIN  VALUE  FOR  THE’ 

PRINT  ’  XTAL  IS  SET  PROPERLY’ 

PRINT  ’IN  LINE  1580’ 

PRINT  ’THE  XTAL  IS  ON  PORT  3’ 

PRINT  ’INC.  POWER  ON  PORT  1’ 

PRINT  ’REF.  POWER  ON  PORT  2* 

INPUT  V* 

IF  V*=CHR*(89)  OR  V*=CHR*(121)  THEN  1460  ELSE  1440 
PRINT  ’ENTER  THE  SCALE  FACTORS  FOR  THE  INC  AND  ’ 
PRINT  ’REF  POWER  METERS- 
INPUT  M(l) 

INPUT  M ( 2 ) 

REM  SETUP  FREQUENCY  SWEEPER 
REM  POWER  LEVEL  RANGE  IN  DBM 
P1=0 
P2  =  1 0 ! 

REM  FREQUENCY  RANGE  IN  GHZ 
F 1  =  .  5 
F2=  1 2 . 5 

REM  NUMBER  OF  FREQUENCY  STEPS  <  514 
N0= 100 

REM  NUMBER  OF  POWER  LEVEL  STEPS  <  10 
N1  =8 

REM  MAX  VOLTAGE  LEVEL  OUT  OF  CRYSTALS 
YMAX* . 1 

REM  GAIN  VALUE  FROM  AMPLIFIER  BOX 
Gl*58 ! 

REM  STEP  SIZES  FOR  LOOPS 
SO*  (F2-F1) / (NO-1) 

Sl=  (P2-P1) / (Nl-1) 

REM  SEND  COMMANDS  TO  SWEEP  OSCILLATOR 
WRT*=’FA’ :CALL  IBWRT (SWEEPER* , WRT*) 

WRT*=STR*(F1) : CALL  IBWRT (SWEEPER% , WRT*) 

WRT*  =  ’  GZ  ’  :  CALL  IBWRT  ( SWEEPER7. ,  WRT* ) 

WRT**  *  FB ’ : CALL  IBWRT ( SWEEPERX , WRT* ) 

WRT*  =  STR*(F2)  :  CALL  IBWRT  (SWEEPER7. ,  WRT*) 

WRT*  *  ’  GZ  ’  :  CALL  I BWRT  ( S  WEEPER7. ,  WRT* ) 

WRT*=  ’  SM’  :  CALL  IBWRT  (SWEEPER7. ,  WRT*) 

WRT*=STR* (FI) : CALL  IBWRT (SWEEPER* , WRT*) 

WRT** ’GZ’ : CALL  IBWRT (SWEEPER* , WRT*) 

WRT»*  ’  SS  ’  :  CALL  IBWRT  ( SWEEPER7. ,  WRT* ) 

WRT**STR«  ( S0»  1 000  (  )  :  CALL  I  BWRT  ( SWEEPER’/. ,  WRT* ) 

WRT*=  "  MZ ’  :  CALL  IBWRT  (SWEEPER*/. ,  WRT*) 

WRT** ’PL’ : CALL  IBWRT  (SWEEPER7. ,  WRT*) 

WRT*=STR* (PI) : CALL  IBWRT (SWEEPERX , WRT*) 

WRT** ’ DM’ : CALL  I BWRT ( SWEEPERX , WRT* ) 

WRT*=’SS’ :CALL  IBWRT (SWEEPERX , WRT*) 

WRT*  =  STR*  (SI)  :  CALL  I  BWRT  ( SWEEPER7. ,  WRT* ) 

WRT**  ’ DM’ :  CALL  IBWRT  (SWEEPER7. ,  WRT* ) 

REM  TAKE  DATA 
CLS 

REM  SET  UP  GRAPH  AXES 
LINE  (25 , 10) - (25 ,210) 

LINE  -(565,210) 


! 

* 
* 


A  A  A  «\ 


MW 


wv 
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MU! 


XTAL(V) 


2  1 0  * ( 1 - ( D ( 3 ) /G 1 ) / YMAX)  ) 
!)*10*-3,D(3)  /G 1 


1825  REM  TURN  ON  RF  SOURCE 

1830  WRT*  =  '  RF  1  '  :  CALL  IBWRT  (SWEEPER7.  ,  WRT*  ) 

1840  FOR  P=P1  TO  P2  STEP  SI 

1850  PRINT  *1.’ ’POWER  LEVEL  (DBM) ’ ’ 

1860  WRITE  #1 ,P 

1870  WRT*  =  '  SM'  :  CALL  I  BWRT  ( SWEEPER’/. ,  WRT* ) 

1880  WRT*=STR* (FI ) :CALL  I BWRT ( SWEEP ER% , WRT* ) 

1890  WRT*  =  ’  GZ  ’  :  CALL  I  BWRT  (SWEEPER*/. ,  WRT* ) 

1900  PRINT  *  1  ,  ’  ’  F ( HZ )  INC(W)  REF ( W)  XTAL(V) 

1910  FOR  F=F 1  TO  F2  STEP  SO 
1920  FOR  1=1  TO  3 
1940  A= I 

1950  Q=M( I ) *5/4095 
2040  REM  GOSUB  2320 
2050  GOSUB  2360 
2060  D ( I ) = Y 

2065  REM  IGNORE  VOLTAGE  LEVELS  FOR  XTAL  AND  REFLECTED  POWER 

2070  IF  1=2  OR  1=3  THEN  2180 

2080  IF  D (I) /M(I) > 1 . 2  THEN  2100 

2090  IF  D(I)/M(I) > .099  THEN  2180 

2100  GOSUB  2710 

2110  GOTO  1940 

2180  NEXT  I 

2185  REM  PLOT  VOLTAGE  VARIATION  OF  XTAL  OUTPUT 

2190  PSET  (( (F-Fl) / (F2-F1) ) *540+25,210* ( 1- (D(3) /Gl) /YMAX) ) 

2200  WRITE  #1 ,F*10'9 ,D ( 1) »10*-3 ,D (2) *10*-3 ,D(3) /G1 

2210  WRT*=’SM"  :CALL  I  BWRT  ( SWEEPER?. ,  WRT* ) 

2220  WRT*  =  *  UP ’ : CALL  IBWRT (SWEEPER* , WRT*) 

2230  NEXT  F 

2240  WRITE  *1,-999,-999,-999,-999 

2260  WRT«= ’PL' :CALL  IBWRT ( SWEEPER* , WRT* ) 

2260  WRT*  =  ’  UP  ’  :  CALL  I  BWRT  ( SWEEPER’/. ,  WRT* ) 

2270  NEXT  P 

2280  WRT*= ’ RFO ’ : CALL  IBWRT (SWEEPER% , WRT* ) 

2290  CALL  I BLOC ( SWEEPER* ) 

2300  CLOSE  *1 
2310  END 

2320  REM  DELAY  SUBROUTINE  *****»***««**»***»***********»**» 1 
2330  FOR  J=1  TO  10 
2340  NEXT  J 
2350  RETURN 

2360  REM  CAMAC  DATA  FETCH  ********************************* 
2370  CALL  IBSIC  (GPIBO*/.) 

2380  WRTS=  ’  CAMAC  ’:  CALL  I BF I ND  ( WRT*  ,  C AMAC’/. ) 

2390  WRT*=CHR* (30) +CHR* (0) +CHR* (17) +CHR* (0) +CHR* (0) +CHR* (0) 
2400  CALL  I BWRT ( CAMAC X , WRT* ) 

24  10  WRT*  =  CHR* (14) +CHR* (0) +CHR* (26) +CHR* (0) +CHR* (0) +CHR* (A) 
2420  CALL  I  BWRT  ( CAMAC  7. ,  WRT*  ) 

2430  WRT*=CHR* (14) +CHR* (0) +CHR* ( 16) +CHR* (0) +CHR* (0) +CHR* (A) 
2440  CALL  I  BWRT  ( CAMAC7. ,  WRT* ) 

2480  WRT*  =  CHR* (14) +CHR* (0) +CHR* (0) 

2490  CALL  I  BWRT  ( CAMAC’/.  ,  WRT*  ) 

2500  CNT%  =  3  :  CALL  I BRDI  (CAMAC7. ,  C27.  ( 0 )  ,  CNT’/.) 

2510  B1’/.  =  C27.  (0) /256 
2  5  20  Y=Q*  (B17,*256+ABS(C27.  (1)  )  +3) 


I 


■/  vv.v. 
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2540 
2550 
2560 
2570 
2600 
2610 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
27  10 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 


WRTS =  CHR$ (14) +  CHRS (0) +  CHRS (24) +CHRS (0) +CHRS (0) +  CHRS (A) 

CALL  IBWRT (CAMAC*/. , WRTS) 

WRTS=CHR$ (30) +CHRS (0) +  CHRS (17) +  CHRS (0) +  CHRS (C) +CHRS (96) 

CALL  IBWRT  (CAMAC*/. ,  WRTS ) 

RETURN 

REM  METER  ZERO  a####****#**#*****####***##*#**##*###*##**## 
WRTS  = " RFO ‘  : CALL  IBWRT ( SWEEPER% , WRTS ) 

LOCATE  1 , 1 

PRINT  'ENTER  A  Y  WHEN  FINISHED  ZEROING  METERS' 

INPUT  V$ 

IF  V*=CHR$ ( 89)  OR  V$=CHRS(121)  THEN  2690  ELSE  2640 
IF  VS  =  CHR$ ( 89 )  OR  V$  =  CHRS(121)  THEN  2700  ELSE  2650 


IF  VS=CHR$ ( 89)  01 
IF  VS  =  CHRS ( 89 )  01 
GOSUB  2320 
RETURN 

REM  SCALE  CHANGE 
IF  1=3  THEN  2800 
K=  I 

LOCATE  1  ,  1 
PRINT  ' 

LOCATE  1  ,  1 
PRINT  'SCALE  SET1 
PRINT  'ENTER  NEW 
INPUT  M ( I ) 

RETURN 


SCALE  SETTING  FOR  PORT' 
ENTER  NEW  VALUE' 


=' ; M( I ) 


v.v.VJ.-.V.V '.'.V 
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DRC . FOR 

C  THIS  PROGRAM  WILL  TAKE  THE  RAW  DATA  OF  A  DIRECTIONAL  COUPLER 
C  CALIBRATION  RUN  AND  PLOT  OUT  THE  ATTENUATION  VS  FREQUENCY. 

C  IF  THIS  PLOT  IS  INCORRECT  THE  INPUT  DATA  FILE  SHOULD  BE  REPEATED. 
DIMENSION  FREQ (514) ,PIN(514) ,POUT(5I4) , ATTN (514) ,TEMP(4) 
DIMENSION  POW (514) 

CHARACTER»70  TEXT  1 , TEXT2 , TEXT3 , TEXT4 , TEXT5 
CHARACTER  FILE*40 , NAME* 10 
CHARACTER)*  20  PTITLE  ,  XLAB  ,  YLAB 
INTEGER  THMI N , THMOUT , DRC 
C  GET  DATA  FILE  OPENED  AND  READY  TO  READ 

WRITE!*,*)  ’ENTER  THE  FILE  NAME  TO  USECIN  SINGLE  QUOTES)’ 
READ!*,*)  FILE 

OPEN (UNIT= 1 , NAME=FILE , STATUS = ’ OLD ’ ) 

READ(1,»)  IFL 

READ ! 1 , * ) TEXT  1 , TEXT5 , 1 PO 

READ ( 1 , » )  TEXT2 

READ ! 1 , * )  DRC 

READ!!,*)  TEXT2 

READ ! 1 , * )  THMIN 

READ ! 1 , * )  TEXT3 

READ ( 1 , * )  THMOUT 

READ ! 1 , * )  TEXT 4 


R£AD(1 , * ) TEXT  1 , TEXT5 , 1 PO 
READ ( 1 , » )  TEXT2 
READ(1,*)  DRC 
READ (I,*)  TEXT2 
READ(1.«)  THMIN 
READ (1 , * )  TEXT3 
READ (1 , * )  THMOUT 
READ ! 1 , * )  TEXT 4 
1  =  1 

READ(1,*)  FREQ! I) .PIN'D  .POUT (I) 

IF  ( (FREQ ! I ) .EQ. -999) .OR. (PIN! I) . EQ . -999) .OR. 
♦GOTO  14 
1  =  1+1 
GOTO  10 
1  =  1-1 

CLOSE (UNIT* 1 ) 

ALL  THE  DATA  IS  READ  IN  NOW  USE  IT 
IF  (I.GT.512)  THEN 

WRITE (*,«)  ’THERE  ARE  TO  MANY  DATA  POINTS’ 
ELSE 

CHECK  FOR  OUTPUT  OF  RAW  DATA  OR  REDUCED  DATA 

WRITE (»,*)  ’ENTER  A  1  FOR  RAW  DATA  OR  A  2  FOR 


( POUT ( I )  . EQ . -999)  ) 


C  CHECK  FOR  OUTPUT  OF  RAW 
WRITE (», *)  ’ENTER  A 
READ!*,*)  K 
IF  (K.EQ.l)  THEN 
C  WE  WILL  CORRECT  OUT  THE 
C  THE  UNCORRECTED  DATA  IS 
C  DATA  IS  RETURNED  IN  TH 
C  CORRECT  PIN 

DO  15  J=  1  , 1 

POW(J) =PIN ( J) 

15  CONTINUE 

CALL  TCOR ( I , FREQ , PO 
DO  16  J=  1  ,  I 

PIN! J) =POW(J) 

16  CONTINUE 
C  CORRECT  POUT 

DO  17  J=  1  ,  I 

POW! J) =POUT ( J) 


FOR  REDUCED  DATA’ 


UT  THE  THERMISTOR  MOUNT  FREQUENCY  RESPONSE 
ATA  IS  FEED  TO  SUBROUTINE  TCOR  CORRECTED 
IN  TH  SAME  ARRAYS. 


POW, THMIN) 


vlv? 


' 


.  •*. 


SS 


Kt 

8 

M 


Es 


23 


17  CONTINUE 

CALL  TCOR ( I , FREQ , POW, THMOUT) 

DO  18  J=  1  , 1 

POUT ( J) =POW(J) 

18  CONTINUE 

C  PLOT  RAW  INCIDENT  AND  PICKOFF  POWER  ON  THE  SAME  SCALE 

C  USE  VAX  GRAPHICS 

C  GRAPHICS  FOR  USE  ON  4010,4014,4025,4105 
WRITE ( *  ,  *  )  ’ENTER  THE  DEVICE  ID#’ 

READ ( * , * )  IDEVIC 

IF  (IDEVIC. EQ. 4014. OR. IDEVIC. EQ. 4105)  THEN 
XMI= . 01 
XMA  = . 9 
YMI= . 01 
YMA  = . 9 

ELSEIF  (IDEVIC. EQ. 4010)  THEN 
XMI= .01 
XMA= .6 
YMI= .01 
YMA= . 6 

ELSEIF  (IDEVIC. EQ. 4025)  THEN 
XMI * . 05 
XMA= .85 
YMI = .05 
YMA= .85 
IDEVI C  =  40 1 0 
ELSE 

WRITE ( * , * )  'DEVICE  ’.IDEVIC,’  NOT  USED’ 

STOP 
END  IF 

CALL  GSTART( IDEVIC) 

CALL  NEWPAG 

CALL  MINMAX(FREQ , I , XWMIN , XWMAX) 

CALL  MINMAX (PIN , I , YMIN , YMAX) 

CALL  MINMAX (POUT, I .YM.YMl) 

IF  (YMIN.LE.YM)  YWMIN=YMIN 
IF  (YMIN.GT.YM)  YWMI N  =  YM 
IF  (YMAX.GE. YM1)  YWMAX= YMAX 
IF  (YMAX.LT. YMi)  YWMAX= YM1 

CALL  MAP (XWMIN, XWMAX, YWMI N.YWMAX, XMI , XMA , YMI , YMA) 
CALL  XLIN (XWMIN, YWMIN, XWMIN, XWMAX, 5) 

CALL  YLIN(XWMIN, YWMIN, YWMIN, YWMAX, 5) 

XLAB= ’ FREQ (HZ) ’ 

YLAB= ’POW(W) ’ 

PTITLE= ’ RAW  DATA’ 

CALL  XLABEL (XLAB) 

CALL  YLABEL ( YLAB) 

CALL  TITLE (PTITLE) 

CALL  CURVE (FREQ, PIN, I , 0 , + 1 ) 

CALL  CURVE (FREQ , POUT , I , 0 , *  1 ) 

CALL  GRSTOP 

C  COMPUTE  ATTENUATION  FOR  UNCORRECTED  DATA  IN  DB 
ELSEIF  (K.EQ.2)  THEN 
C  REMOVE  THEMSTOR  CAL  FACTORS 
DO  21  J= 1 , I 
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POW(J) =PIN(J) 

21  CONTINUE 

CALL  TCOR ( I , FREQ , POW, THMIN) 

DO  22  J=  1  , 1 

PIN ( J) =POW( J) 

22  CONTINUE 

DO  23  J=  1  , 1 

POW(J) =POUT (J) 

23  CONTINUE 

CALL  TCOR ( I , FREQ , POUT , THMOUT) 

DO  24  J=1  , 1 

v  POUT ( J) =POW( J) 

■<.  24  CONTINUE 

DO  30  J= 1 , I 

ATTN ( J) = 10 . #AL0G1G (PIN ( J) /POUT (J) ) 
30  CONTINUE 
C  PLOT  ATTN  VS  FREQ. 

-  C  USE  VAX  GRAPHICS 

■.  C  GRAPHICS  FOR  USE  ON  4010,4014,4025,4105 

WRITE ( * , * )  ’ENTER  THE  DEVICE  ID«* 
j  READ ( * , «)  IDEVIC 

IF  (IDEVIC. EQ. 4014. OR. IDEVIC. EQ. 4105) 


THEN 


(IDEVIC. EQ. 4025)  THEN 


XMI =.01 
XMA= . 9 
YMI =.01 
YMA= .9 

ELSEIF  (IDEVIC. EQ. 4010)  THEN 
XMI =.01 
XMA= .6 
YMI  =  .  0 1 
YMA=.6 

ELSEIF  (IDEVIC. EQ. 4025)  THEN 
XMI = . 05 
XMA= .85 
YMI = .05 
YMA= .85 
IDEVIC=40 10 
ELSE 

WRITE ( * , * )  ’DEVICE  ', IDEVIC,’  NOT  USED’ 

STOP 
END  IF 

CALL  GSTART( IDEVIC) 

CALL  NEWPAG 

CALL  MINMAX (FREQ , I , XWMIN , XWMAX) 

CALL  MINMAX (ATTN, I , YWMIN , YWMAX) 

CALL  MAP (XWMIN, XWMAX, YWMIN, YWMAX, XMI , XMA , YMI , YMA) 
CALL  XLIN (XWMIN, YWMIN, XWMIN, XWMAX, 5) 

CALL  YLIN (XWMIN, YWMIN, YWMIN, YWMAX, 5) 

XLAB= ’ FREQ (HZ) ’ 

YLAB= ’ ATTN (DB) ’ 

PTITLE= ’ FREQ .  RESPONSE’ 

CALL  XLABEL (XLAB) 

CALL  YLABEL (YLAB) 

CALL  TITLE (PTITLE) 

CALL  CURVE (FREQ, ATTN, I ,0,+l) 


s.  ■  _  v  y. 


l,'*  ***.>,fc Jl* 


COUPLER  ID*= 


FREQ (HZ) ’ ,20X, ’ATTN(DB) 
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CALL  GRSTOP 

C  CREATE  FILE  FOR  LOOKING  AT  ATTN  VALUES 

OPEN (UNIT= 1 ,  F I LE= ’ DRATTN . PAT  *  , STATUS= ’ NEW' ) 

WRITE ( 1 , * )  ’’’FLAG  SET  TO  1  MEANS  DATA  CORRECTED’ 
WRITE ( 1 , » )  1 

WRITE ( 1 , * )  ’’’DIR.  COUPLER  ID*  =  ’” 

WRITE ( 1 , *  >  DRC 
WRITE ( 1 , « )  I 
WRITE(1 ,40) 

40  FORMAT ( ’  ’,’’’ FREQ (HZ) ’, 20X, ’ ATTN (DB) ’’’ ) 

DO  50  J= 1 , I 

WRITE (1,45)  FREQ ( J) , ATTN ( J ) 

45  FORMAT ( ’  ’  , 1 PE  1 1 . 3 , 1 7X , 1  PE  1 1 . 3 ) 

50  CONTINUE 

CLOSE (UNIT= 1 ) 

END  IF 
END  IF 
STOP 
END 
C 

C  THIS  SUBROUTINE  TCOR 

SUBROUTINE  TCOR (J , FREQ , POW , ITHMIN) 

DIMENSION  FREQ ( J) , POW{ J) 

REAL « 8  NU ( 50 ) ,CFAC(50) ,B(50) ,C(50) ,D(50) , FRQ , CAL 
CHARACTER  FNAME«20 

C  GET  THE  DATA  FILE  THERMISTOR  MOUNT  ITHMIN 
IF  (ITHMIN. EQ. 21593)  FNAME= ’ 21593 . MNT ’ 


21594 .MNT’ 
18832. MNT’ 
21627. MNT’ 
51009. MNT’ 
50918. MNT’ 
19131 .MNT’ 


IF  (ITHMIN. EQ. 21594)  FNAME= ' 2 1 594 . MNT ’ 

IF  (ITHMIN. EQ. 18832)  FNAME* ' 18832 . MNT ’ 

IF  (ITHMIN. EQ. 21627)  FNAME= ' 2 1 627 . MNT ’ 

IF  (ITHMIN. EQ. 51009)  FNAME= ’ 51009 . MNT ’ 

IF  (ITHMIN. EQ. 50918)  FNAME= ' 509 18 . MNT ’ 

IF  (ITHMIN. EQ. 19131)  FNAME= ’ 19 13 1 . MNT ’ 

OPEN (UNIT=2 , FILE=FNAME , STATUS = ’OLD’ ) 

GET  THE  NUMBER  OF  DATA  PAIR 
READ ( 2 , * )  MAX 
DO  10  1=1, MAX 

READ (2,2)  NU ( I ) , CFAC ( I ) 

FORMAT (2D  10. 2) 

)  CONTINUE 

CLOSE ( UNI T=  2 ) 

GET  THE  INTERPOLATION  ARRAYS  FILLED 

CALL  SPCOEF ( MAX , NU , CFAC , B , C , D , IFLAG) 

IF  (IFLAG.NE.O)  THEN 

WRITE ( » , » )  ’ERROR  IN  INTERPOLATON ,  ’, IFLAG 
END  IF 

CORRECT  ARRAY  POW 

ARRAY  NU  IS  IN  GHZ,  AND  ARRAY  FREQ  IS  IN  HZ  SO  WE  MUST  TAKE 

HE  DIFFERENT  UNITS 
DO  20  1=1 , J 

FRQ=DBLE(FREQ(I) ) * ( 1 .D-9) 

CALL  S VALUE (MAX, NU , CFAC , B , C , D , FRQ , CAL , IFLAG) 

IF  ( IFLrtG. NE . 0)  THEN 

WRITE (»,»)  ’ERROR  IN  EVALUATION,  ’  , I  FLAG , FREQ ( I ) 
ENDIF 


IFLAG 


vy/.vNVv: 
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I 


v>>.* 
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POW(I) = (POW(I) ) /CAL 
20  CONTINUE 

C  THE  ARRAY  POW  IS  NOW  CORRECTED  FOR  THE  THERMISTOR  MOUNT  USED 
RETURN 
END 


27 
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PROGRAM  DB 

C  THIS  PROGRAM  WILL  READ  A  FILE  GENERATED  BY  THE  ATTN  PROGRAM  ON  THE 
C  4051  CALIBRATION  SYSTEM 

DIMENSION  FREQ (514) ,P0W(514,3) , ITHERM(3) ,IDRC(2) ,TEMP(514) 
CHARACTER  LABEL  1*20, LABEL2#50 , XLAB« 1 5 . YLAB# 1 5 , PTITLE* 1 5 
CHARACTER  FNAME«40 
C  GET  THE  DATA  FILE  NAME  TO  LOOK  AT 

WRITE (*,*)  ’ENTER  THE  DATA  FILE  NAME’ 

READ ( * , * )  FNAME 

OPEN (UNIT= 1 ,FILE=FNAME, STATUS =’ OLD’ ) 

READ(1,*)  LABEL  1 
READ  (1,0  IFLAG 
READ ( 1 , * )  LABEL  1 
READ ( 1 , * )  ISMPLE 
READ  (1,0  LABEL  1 
READ  (1,0  IDRC(l) 

READ ( 1 , * )  LABEL  1 
READ  (1,0  IDRC  ( 2 ) 

READ  (1,0  LABEL  1 
READ  (1,0  I  THERM  ( 1 ) 

READ  (1,0  LABEL  1 
READ  (1,0  ITHERM(2) 

READ  (1,0  LABEL  1 
READ  (1,0  I  THERMO) 

READd.O  LABEL2 
1  =  1 

5  READd.O  FREQ  ( I )  ,  POW  (1,1)  ,  POW(  1,2),  POW  (1,3) 

IF (FREQ (I)  . EQ. -999. OR. POW (I , 1)  . EQ . -999 . OR . POW( I ,2)  .EQ.-999.0R. 
+POW( I ,3) . EQ . -999)  GOTO  10 
1  =  1  +  1 
GOTO  5 

10  1=1-1 

CLOSE (UNIT= 1 ) 

IF  (I.GT.512)  THEN 

WRITE (»,*)  ’TO  MANY  DATA  POINTS’ 

STOP 
END  IF 

C  REMOVE  THERMISTOR  RESPONSE  FROM  POW(I,J)  FOR  J=l,2,3 
DO  15  K=1 ,3 
DO  12  L=  1  ,  I 

TEMP (L) =POW(L,K) 

12  CONTINUE 

ITH=ITHERM(K) 

CALL  TCOR ( I , FREQ , TEMP , ITH) 

DO  14  L=  1  .  I 

POW ( L , K) =TEMP  (L) 

14  CONTINUE 

15  CONTINUE 

C  REMOVE  THE  DIRECTIONAL  COUPLER  RESPONSE  FROM  POW(L,K)  FOR  K=l,2 
DO  25  K=  1  ,  2 
DO  20  L= 1 , I 

TEMP (L) =POW(L , K) 


i 

& 
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20  CONTINUE 

ID= IDRC (K) 

CALL  DCCOR ( I , FREQ , TEMP , ID) 

DO  24  L=  1  , 1 

P0W(L , K) =TEMP (L) 

24  CONTINUE 

25  CONTINUE 

C  NOW  THE  POWER  VALUES  HAVE  BEEN  CORRECTED  FOR  THE  FREQUENCY  RESPONSE 

C  OF  THE  THERMISTOR  MOUNTS  AND  THE  DIRECTIONAL  COUPLERS  USED  TO  TAKE 

C  THE  VARIOUS  POWER  VALUES. 

C  COMP  ATTENUATION,  ATTN= (POW (L , 1) -P0W(L,2) ) /POW ( L , 3 ) 

29  WRITE (*,*)  ’ENTER  A  1  IF  YOU  WANT  THE  ATTENUATION  COMPUTED’ 
WRITE!*, »)  ’WITH  ONLY  THE  INCIDENT  POWER’ 

WRITE (»,*)  ’ENTER  A  0  IF  YOU  WANT  THE  ATTENUATION  COMPUTED’ 
WRITE!*,*)  ’USING  THE  NET  POWER’ 

READ!*,*)  IOPT 

IF  (IOPT. NE. 1 .AND. IOPT. NE.O)  GOTO  29 
DO  30  L=  1  ,  I 

IF  (IOPT.EQ.O)  THEN 

ATTN= (POW(L , 1) -POW(L ,2) ) /POW(L ,3) 

C  TEMP (L) =0 .  MEANS  THEIR  IS  AN  ERROR  IN  THE  DATA 
IF  (ATTN.LE.l.)  THEN 
TEMP (L) =0. 

ELSE 

TEMP (L) = 10 . *AL0G10 ( ATTN) 

END  IF 
ELSE 

ATTN=P0W(L, 1) /P0W(L,3) 

IF  (ATTN.LE.l)  THEN 
TEMP !L) =0. 

ELSE 

TEMP (L) = 10. *AL0G10 (ATTN) 

END  IF 
ENDIF 

30  CONTINUE 

C  MAKE  A  GRAPH  OF  THE  FREQUENCY  RESPONSE 
XLAB= ’FREQ (HZ) ’ 

YLAB= ’ ATTN(DB) ’ 

C  USE  VAX  GRAPHICS 

C  GRAPHICS  FOR  USE  ON  4010,4014,4025,4105 
WRITE ( * , * )  ’ENTER  THE  DEVICE  ID*’ 

READ ( * , * )  IDEVIC 

IF  (IDEVIC. EQ. 4014. OR. IDEVIC. EQ. 4105)  THEN 
XMI= .01 
XMA= .9 
YMI =.01 
YMA= .9 

ELSEIF  (IDEVIC. EQ. 4010)  THEN 
XMI= .01 
XMA= . 6 
YMI = .01 
YMA= . 6 

ELSEIF  (IDEVIC . EQ . 4025)  THEN 
XMI = . 05 
XMA= . 85 
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YMA  = .85 
IDEVIC=40 1 0 
ELSE 

WRITE ( *  ,  * )  ’DEVICE  ’.IDEVIC,’  NOT  USED’ 

STOP 
END  IF 

CALL  GSTART( IDEVIC) 

CALL  NEWPAG 

CALL  MINMAX (FREQ , I , XWMIN , XWMAX) 

CALL  MINMAX (TEMP , I , YWMI N , YWMAX) 

CALL  MAP (XWMIN, XWMAX, YWMI N, YWMAX , XMI , XMA , YMI , YMA) 
CALL  XLIN (XWMIN, YWMIN, XWMIN .XWMAX, 5) 

CALL  YLIN (XWMIN , YWMIN , YWMIN , YWMAX ,5) 

IF  (IOPT.EQ.O)  PTITLE= ’ NET  POWER’ 

IF  (IOPT.EQ.l)  PTITLE= ’ INC .  POWER’ 

CALL  XLABEL (XLAB) 

CALL  YLABEL ( YLAB) 

CALL  TITLE (PTITLE) 

CALL  CURVE (FREQ, TEMP, I ,0,+I) 

CALL  GRSTOP 

C  COMPUTE  ATTENUATION  FOR  UNCORRECTED  DATA  IN  DB 
C  CREATE  FILE  FOR  SAVING  DATA 

OPEN (UNIT= 1 ,FILE=' SAMPLE. DAT’ , STATUS* ’ NEW’ ) 


WRITE ( I , * )  ’’’FLAG  SET  TO  1  MEANS  DATA  IS  CORRECTED ' 
WRITE ( 1 , « )  1 

WRITE (1,*)  ’’’SAMPLE  ID  TAG’” 

WRITE ( 1 , * )  ISMPLE 
WRITE ( 1 , * )  I 
WRITE (1 ,40) 

>  FORMAT ( ’  ’  ,  ”  ’FREQ(HZ)  ’  , 19X, ’ATTN(DB)  ”  ’) 

DO  50  L= 1  , 1 

WRITE (1,45)  FREQ (L) .TEMP (L) 

>  FORMAT  (  '  ’  ,  1PEI 1 .3,  I7X,  IPEU  .3) 

)  CONTINUE 

CLOSE (UNIT=1) 

STOP 
END 

SUBROUTINE  TCOR 

SUBROUTINE  TCOR (J , FREQ , TEMP , ITHERM) 

DIMENSION  FREQ (J) , TEMP (J) 

REAL*8  NU ( 50 )  ,CFAC(50)  ,B(50)  ,C(50)  ,D(50)  ,  FRQ , CAL 
CHARACTER»20  FNAME 


FREQ (HZ)  ’  , 19X, ' ATTN (DB)  ’  ’  ’ ) 

FREQ (L) , TEMP (L) 

1PEI 1 .3,  I7X,  IPEU  .3) 


IF  (ITHERM. EQ. 21593)  FNAME= ’ 2 1 593  .  MN 
IF  (ITHERM. EQ. 21594)  FNAME= ’ 2 1 594 . MN 
IF  (ITHERM. EQ. 18832)  FNAME = ’ 1 8832 . MN 
IF  (ITHERM. EQ. 21627)  FNAME = ’ 2 1 627 . MN 
IF  (ITHERM. EQ. 51009)  FNAME= ’ 5 1 009  .  MN 
IF  (ITHERM. EQ. 50918)  FNAME= ’ 509 1 8 . MN 
IF  (ITHERM. EQ. 19131)  FNAME= ’ 1 9 1 3 1 . MN 
OPEN ( UNIT=  2 , FI LE  =  FNAME .STATUS = ’OLD’ ) 
GET  THE  NUMBER  OF  DATA  POINTS 
READ ( 2 , # )  MAX 
DO  10  1=1, MAX 


MNT  ’ 
21594 . MNT’ 
18832 . MNT’ 
21627 .MNT’ 
51009. MNT’ 
50918 . MNT’ 
19131. MNT ’ 
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READ (2,2)  NU ( I ) ,CFAC(I) 

2  FORMAT (2D  10. 2) 

10  CONTINUE 

C  CORRECT  ARRAY  TEMP  BY  INTERPOLATING  ARRAY  CFAC 
CALL  SPCOEF ( MAX , NU , CFAC . B , C , D , I  FLAG) 

IF  ( I FLAG . NE . 0)  THEN 

WRITE (*,*) ’ERROR  IN  INTERPOLATION  ’ , I FLAG 
END  IF 

DO  20  1=1  ,  J 

FRQ=DBLE (FREQ ( I ) ) * (1 .D-9) 

CALL  SVALUE( MAX, NU, CFAC ,  B , C , D , FRQ , CAL , I FLAG) 

IF  (IFLAG.NE.O)  THEN 

WRITE (#,*) 'ERROR  IN  EVALUATION  ’ , IFLAG , FREQ ( I ) 

ENDIF 

TEMP (I) =TEMP ( I ) /CAL 
20  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  DCCOR(J , FREQ , TEMP , ID) 

DIMENSION  FREQ (J) , TEMP (J) 

REAL #8  NU (514) ,PCK0FF(514) ,B(514) ,C(514) ,D(514) , FRQ , VAL 
CHARACTER  FNAME»20 , LABEL»50 
IF  (ID.EQ.l)  FNAME= ’ DC  1 . DAT ’ 

IF  (ID.EQ.2)  FNAME= ’ DC2 . DAT ' 

IF  (ID.EQ.3)  FNAME= ’ DC3 . DAT ’ 

IF  (ID.EQ.4)  FNAME= ’ DC4 . DAT ’ 

OPEN (UNIT=2 , FILE=FNAME , STATUS* ' OLD ’ ) 

READ ( 2 , » )  LABEL 
READ (2,*)  I 
READ ( 2 , * )  LABEL 
READ ( 2 , * )  I 
READ ( 2  ,  » )  MAX 
READ  (2  ,  « )  LABEL 
DO  5  1=1, MAX 

READ ( 2 , » )  A.Z 
NU ( I ) =DBLE (A) 

PCKOFF ( I ) =DBLE (Z) 

5  CONTINUE 

C  INTERPOLATE  THE  DATA 

CALL  SPCOEF (MAX , NU , PCKOFF , B , C , D , IFLAG) 

IF  (IFLAG.NE.O)  THEN 

WRITE ( * , * )  ’ERROR  IN  INTERPOLATION’ , IFLAG 
ENDIF 

C  REMOVE  DIR.  COUPLERSE 
DO  10  1=1 , J 

FRQ=DBLE (FREQ ( I ) ) 

CALL  S VALUE ( MAX , NU , PCKOFF , B , C , D , FRQ , VAL , I  FLAG) 

IF  (IFLAG.NE.O)  THEN 

WRITE (»,«) 'ERROR  IN  PICKOFF  EVALUATION’ , IFLAG , FREQ ( I ) 
ENDIF 

TEMP ( I ) =VAL»TEMP ( I ) 

10  CONTINUE 
RETURN 
END 
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XTAL . FOR 

PROGRAM  XTAL 

C  THIS  PROGRAM  WILL  READ  A  FILE  GENERATED  BY  THE  XTAL  PROGRAM  ON  THE 
C  4051  CALIBRATION  SYSTEM 
C  FREQ(I.J)  HOLDS  THE  FREQUENCY 
C  POWIN(I.J)  HOLDS  THE  INCIENT  PICKOFF  POWER 
C  POWREF ( I , J )  HOLDS  THE  REFLECTED  PICKOFF  POWER 
C  VOLT ( I , J)  HOLDS  THE  XTAL  VOLTAGE 
C  I  IS  THE  COUNTER  USED  FOR  FREQUENCY 
C  J  IS  THE  COUNTER  USED  FOR  POWER  APPLIED 
PARAMETER (IMAX=5 14 , IPMAX= 10) 

DIMENSION  FREQ ( I MAX , IPMAX)  , POWI N ( I  MAX , I PMAX)  . POWREF ( I  MAX . I PMAX) 
DIMENSION  VOLT (I MAX, IPMAX) ,IDRC(2) , ITHERM ( 2 ) .PWR(IPMAX) 
DIMENSION  TEMPP(IMAX) .TEMPF(IMAX) 

CHARACTER  LABEL  1 *20 , LABEL2» 50 , XLAB# 1 5 , YLAB» 1 5 .PTITLE*  15 
CHARACTER  FNAME#40 
C  GET  THE  DATA  FILE  NAME  TO  LOOK  AT 

WRITE ( * , « )  ’ENTER  THE  DATA  FILE  NAME’ 

READ ( * , » )  FNAME 

OPEN (UNIT= 1 , FILE=FNAME , STATUS = ’ OLD ’ ) 

READ ( 1 , * )  LABEL  1 
READ ( 1 , » )  I FLAG 
READ ( 1 , * )  LABEL  1 
READ(1,*)  I XTAL 
READ ( 1 , « )  LABEL  1 
READ ( 1 , * )  IDRC(l) 

READ ( 1 , » )  LABEL  1 
READ ( 1 , »)  IDRC (2) 

READ ( 1 , * )  LABEL  1 
READ(1,*)  ITHERM ( 1 ) 

READ ( 1 , « )  LABEL  1 
READd,*)  ITHERM(2) 

J  =  0 

1  J  =  J+1 

READd  ,*  ,END=10)  LABEL  1 
C  PWR  READS  THE  APPLIED  POWER  LEVEL 
READ (1  ,  » )  PWR(J) 

READ(1,»)  LABEL2 
1  =  1 

5  READ(1, *)  FREQ ( I , J) , POWIN ( I , J) , POWREF (I , J) , VOLT ( I , J) 

IF (FREQ (I ,J) .EQ. -999. OR. POWIN (I ,J) . EQ . -999 . OR . POWREF (I ,J) . EQ . 
+-999  .OR. VOLTCI , J) . EQ . -999)  GOTO  1 
1  =  1  +  1 
GOTO  5 

10  CLOSE (UNIT* 1 ) 

1  =  1-1 
J  =  J-  1 

IF  (I.GT.512)  THEN 

WRITE ( * , * )  ’TO  MANY  DATA  POINTS’ 

STOP 
END  IF 

C  REMOVE  THERMISTOR  RESPONSE  FROM  POWIN(I.J) 

C  AND  POWREF ( I, J)  FOR  J=l,2 
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DO  15  K=  1  ,  J 
DO  12  L=  1  ,  I 

TEMPF (L) =  FREQ (L , K) 

TEMPP (L) =POWI N  ( L  ,  K) 

12  CONTINUE 

ITH=ITHERM( 1) 

CALL  TCOR (I .TEMPF, TEMPP . ITH) 

DO  14  L= 1 , I 

POWI N ( L , K) =TEMPP (L) 

14  CONTINUE 

15  CONTINUE 

DO  19  K=  1  ,  J 

DO  17  L= 1 , I 

TEMPF (L) =FREQ(L,K) 

TEMPP (L) =POWREF (L , K) 

17  CONTINUE 
ITH= ITHERM( 2 ) 

CALL  TCOR ( I , TEMPF .TEMPP , ITH) 

DO  18  L= 1  , 1 

POWREF (L , K) =TEMPP (L) 

18  CONTINUE 

19  CONTINUE 

C  REMOVE  THE  DIRECTIONAL  COUPLER  RESPONSE 
DO  28  K=  1  ,  J 
DO  20  L=  1  ,  I 

TEMPF (L)=FREQ(L,K) 

TEMPP (L) =POWIN (L ,K) 

20  CONTINUE 
ID= IDRC ( 1 ) 

CALL  DCCOR( I .TEMPF, TEMPP. ID) 

DO  24  L=  1  , 1 

POWIN(L.K) =TEMPP (L) 

24  CONTINUE 

DO  25  L=  1  , 1 

TEMPF (L) =  FREQ (L , K) 

TEMPP (L) =POWREF(L,X) 

25  CONTINUE 
ID=IDRC (2) 

CALL  DCCOR(I .TEMPF, TEMPP , ID) 

DO  26  L= 1  , 1 

POWREF (L.K) =TEMPP (L) 

26  CONTINUE 

28  CONTINUE 
C 

29  WRITE (»,*)  ’ENTER  A  1  IF  YOU  WANT  NET  POWER  COMPUTED’ 

WRITE  (.*,*)  ’ENTER  A  0  IF  YOU  WANT  INCIDENT  POWER  ONLY’ 

READ ( » , * )  IOPT 

IF  ( IOPT. NE. 1 .AND. IOPT. NE.O)  GOTO  29 
C  NOW  THE  POWER  VALUES  HAVE  BEEN  CORRECTED  FOR  THE  FREQUENCY  RESPONSE 
C  OF  THE  THERMISTOR  MOUNTS  AND  THE  DIRECTIONAL  COUPLERS  USED  TO  TAKE 
C  THE  VARIOUS  POWER  VALUES. 

C  THE  NET  POWER  WILL  BE  STORED  IN  POWIN,  AND  DATA  FILES  FOR  BOTH  THE 
C  APPLIED  POWER  AND  THE  INDUCED  XTAL  RESPONSE  WILL  BE  GENERATED. 

C  SAVE  NET  POWER 

IF  (IOPT.EQ.l)  THEN 
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DO  40  K=  1  ,  J 
DO  30  L= 1 , I 

POWIN(L.K) =  POWI N ( L , K) -POWREF(L.K) 
30  CONTINUE 

40  CONTINUE 

END  IF 


C 

C  SET  UP  FOR  GRAPHICS 

C  USE  VAX  GRAPHICS 

C  GRAPHICS  FOR  USE  ON  4010,4014,4025,4105 

WRITE (*,*)  ’ENTER  THE  DEVICE  ID* ’ 

READ ( * , * )  IDEVIC 

IF  (IDEVIC . EQ . 4014 . OR. IDEVIC . EQ . 4 105)  THEN 
XMI = .01 
XMA= . 9 
YMI= . 01 
YMA= .9 

ELSEIF  (IDEVIC. EQ. 4010)  THEN 
XMI = .01 
XMA= .6 
YMI =.01 
YMA= .6 

ELSEIF  (IDEVIC .EQ. 4025)  THEN 
XMI = .05 


XMA= .85 
YMI = .05 
YMA= .85 
IDEVIC  =  40  10 
ELSE 

WRITE (  »  ,  » )  'DEVICE  ’.IDEVIC,’  NOT  USED’ 

STOP 
END  IF 

C  PLOT  POWER  VS  FREQ. 

CALL  GSTART( IDEVIC) 

CALL  NEWPAG 
XWMIN=FREQ ( 1 , 1) 

XWMAX=  XWMI N 
YWMIN=POWIN( 1 . 1) 

YWMAX=YWMIN 
DO  70  K= 1 , J 
DO  65  L= 1 , I 

TEMPF(L) =FREQ (L , K) 

TEMPP(L) =POWIN (L , K) 

65  CONTINUE 

CALL  MINMAX (TEMPF , I , XI , X2 ) 

CALL  MINMAX (TEMPP , I , Y1 , Y2 ) 

IF  (XWMIN.GT.X1)  XWMI N= XI 
IF  (XWMAX.LT.X2)  XWMAX=X2 
IF  (YWMIN.GT. Yl)  YWMIN=Y1 
IF  (YWMAX.LT. Y2)  YWMAX=Y2 
70  CONTINUE 

CALL  MAP ( XWMI N , XWMAX , YWMI N , YWMAX , XMI , XMA , YMI , YMA) 
CALL  XL IN (XWMIN , YWMIN , XWMIN , XWMAX , 5 ) 

CALL  YLIN (XWMIN, YWMIN, YWMIN, YWMAX, 5) 

XLAB= ’ FREQ. (HZ) ’ 


34 


% 


V!  ' 


-  -  V." 


*.■**.’ 


AFWL-TR-87-119 


YLAB= ' POWER  (W)  ’ 

IF  (I0PT.EQ.1)  PTITLE= ’ NET  POWER  RESPONSE’ 

IF  (IOPT.EQ.O)  PTITLE= ’ INCIDENT  POWER  RESPONSE’ 
CALL  XLABEL ( XLAB) 

CALL  YLABEL (YLAB) 

CALL  TITLE (PTITLE) 

DO  76  K= 1 , J 


DO  74  L=  1  ,  I 

TEMPF(L) =FREQ (L , K) 

TEMPP(L) =  POWI N ( L , K ) 

74  CONTINUE 

CALL  CURVE (TEMPF .TEMPP , I ,0 ,0) 

76  CONTINUE 

CALL  GRSTOP 
78  READ ( * , * )  M 

IF  (M.NE. 1)  GOTO  78 

C  PLOT  VOLTAGE  RESPONSE  TO  APPLIED  POWER 
CALL  GSTART ( IDEVIC) 

CALL  NEWPAG 
XWMIN=FREQ ( 1 , 1) 

XWMAX=XWMIN 
YWMIN=VOLT ( 1 , 1) 

YWMAX= YWMIN 
DO  85  K=  1  ,  J 
DO  80  L=  1  ,  I 

TEMPF (L) =FREQ (L  ,  K) 

TEMPP(L) = VOLT (L , K) 

80  CONTINUE 

CALL  MI NMAX( TEMPF , I ,X1 ,X2) 

CALL  MI NMAX ( TEMPP ,1 , Y 1 , Y2 ) 

IF  (XWMIN.GT.X1)  XWMIN=X1 
IF  (XWMAX.LT. X2)  XWMAX=X2 
IF  (YWMIN. GT.Y1)  YWMIN=Y1 
IF  (YWMAX.LT. Y2)  YWMAX=Y2 
85  CONTINUE 

CALL  MAP (XWMIN.XWMAX, YWMIN, YWMAX.XMI , XMA , YMI , YMA) 
CALL  XL IN (XWMIN , YWMIN , XWMIN , XWMAX , 5) 

CALL  YLIN (XWMIN, YWMIN, YWMIN, YWMAX, 5) 

XLAB= 'FREQ. (HZ) ’ 
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YLAB= ’VOLTS’ 

PTITLE= ’ XTAL  FREQ. 

CALL  XLABEL (XLAB) 

CALL  YLABEL (YLAB) 

CALL  TITLE (PTITLE) 

DO  90  K= 1 , J 
DO  88  L=  1  ,  I 

TEMPF (L) =FREQ (L , K) 
TEMPP (L) = VOLT (L , K) 
CONTINUE 

CALL  CURVE (TEMPF .TEMPP 
CONTINUE 
CALL  GRSTOP 


1,0,0) 


C 

C  CREATE  FILE  FOR  SAVING  DATA 

OPEN ( UNIT= 1 , FI LE= ’XTALP.DAT’ , STATUS = ’NEW’ ) 
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0PEN(UNIT=2,FILE=’ XTALV . DAT ' . STATUS = ’ NEW’ ) 

WRITE (1,»)  ’  ’  ’ XTAL  ID  TAG ’  ’  * 

WRITE (2,*)  ’’’XTAL  ID  TAG’’’ 

WRITE ( 1 , * )  I XTAL 
WRITE (2,*)  I XTAL 

WRITE (1,«)  ’’’SWEEPER  POWER  OUTPUT’ ’’ 

WRITE (2,»)  ’’’XTAL  OUTPUT  VOLTAGE’’’ 

WRITE ( 1 , * )  I.J 
WRITE  (2  ,  * )  I.J 
WRITE (1 ,91)  (PWR(K)  ,K=I  .  J) 

WR I T  E ( 2 ,91)  (PWR(K)  ,K=1  , J) 

91  FORMAT ( ’  ’ , ’ ’ ‘ FREQXPOW’ ’ ’ ,2X, 10 (F8 . 2 ,5X) ) 

DO  95  L=  1  ,  I 

WRITE (2. 92)  FREQ ( L  ,  1)  .  (V0LT(L,K)  ,K=1  ,J) 

WRITE (1 .92)  FREQ ( L ,  1)  ,  (POWIN(L.K)  ,K=1  ,J) 

95  CONTINUE 

92  FORMAT ( ’  ' , 1 1 ( 1PE1 1 . 3 , 2X) ) 

CLOSE ( UNIT= 1 ) 

CLOSE (UNIT=2) 

STOP 
END 

SUBROUTINE  TCOR 

SUBROUTINE  TCOR ( J , TEMPF , TEMPP , ITHERM) 

DIMENSION  TEMPF (J) .TEMPP (J) 

REAL *8  NU ( 50 ) ,CFAC(50) ,B(50) ,C(50) ,D(50) .FRQ.CAL 
CHARACTER* 20  FNAME 

IF  (ITHERM. EQ. 21593)  FNAME= ’ 2 1593 . MNT ’ 

IF  (ITHERM. EQ. 21594)  FNAME= ’ 2 1 594 . MNT ’ 

IF  (ITHERM. EQ. 18832)  FNAME= ’ 18832 . MNT ’ 

IF  (ITHERM. EQ. 21627)  FNAME= ’ 2 1627 . MNT ’ 

IF  (ITHERM. EQ. 51009)  FNAME= ’ 5 1 009 . MNT ’ 

IF  (ITHERM. EQ. 50918)  FNAME= ’ 509 18 . MNT ’ 

IF  (ITHERM. EQ. 19131)  FNAME= ’ 19131 . MNT ’ 

OPEN (UNIT=2 , FILE=FNAME , STATUS = ’ OLD ' ) 

C  GET  THE  NUMBER  OF  DATA  POINTS 
READ (2,*)  MAX 
DO  10  1=1, MAX 

READ (2, 2)  NU ( I ) ,CFAC(I) 

2  FORMAT (2D10 . 2) 

10  CONTINUE 

CLOSE (UNIT=2) 

C  CORRECT  ARRAY  TEMPP  BY  INTERPOLATING  ARRAY  CFAC 
CALL  SPCOEF ( MAX , NU , CFAC , B , C , D , I FLAG) 

IF  ( I  FLAG . NE . 0 )  THEN 

WRITE (*,*) ’ERROR  IN  INTERPOLATION  ’ , IFLAG 
END  IF 

DO  20  1  =  1  .J 

FRQ=DBLE (TEMPF ( I ) ) * ( 1 . D- 9 ) 

CALL  SVALUE ( MAX , NU , CFAC , B , C , D , FRQ , CAL , I  FLAG) 

IF  (IFLAG. NE.O)  THEN 

WRITE (»,*)  ’ERROR  IN  EVALUATION  ’  . I  FLAG , TEMPF ( I ) 
END  I  F 

TEMPP (I) =TEMPP(I) /CAL 
20  CONTINUE 
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RETURN 

END 


SUBROUTINE  DCCOR( J .TEMPF .TEMPP , ID) 

DIMENSION  TEMPF (J) , TEMPP (J) 

REAL #8  NU (514) ,PCKOFF(514) ,B(5I4) ,C(514) ,D(514) , FRQ , VAL 
CHARACTER  FNAME*20 ,LABEL*50 
IF  (ID.EQ.l)  FNAME= ’ DC  1 .DAT’ 

IF  (ID.EQ.2)  FNAME= ’ DC2 . DAT ’ 

IF  (ID.EQ.3)  FNAME= ’ DC3 . DAT ’ 

IF  (ID.EQ.4)  FNAME= ’ DC4 . DAT ' 

OPEN (UNIT=2 , FILE=FNAME , STATUS =’ OLD ’ ) 

READ (2,#)  LABEL 
READ (2,*)  I 
READ ( 2 , » )  LABEL 
READ ( 2 , * )  I 
READ (  2  ,  «  )  MAX 
READ ( 2 , * )  LABEL 
DO  5  1=1, MAX 

READ (2,*)  A , Z 
NU ( I ) =DBLE (A) 

PCKOFF ( I ) =DBLE (Z) 

5  CONTINUE 

CLOSE (UNIT=2) 

C  INTERPOLATE  THE  DATA 

CALL  SPCOEF (MAX , NU , PCKOFF , B . C , D , I FLAG) 

IF  ( I FLAG. NE . 0)  THEN 

WRITE ( * , » )  ’ERROR  IN  INTERPOLATION*  , IFLAG 
END  IF 

C  REMOVE  DIR.  COUPLERSE 
DO  10  1=1, J 

FRQ  =  DBLE (TEMPF ( I )  ) 

CALL  SVALUE ( MAX. NU, PCKOFF, B,C ,D, FRQ, VAL, IFLAG) 

IF  (IFLAG. NE.O)  THEN 

WRITE («,*) 'ERROR  IN  PICKOFF  EVALUATION' , IFLAG , TEMPF ( I ) 
END  IF 

TEMPP ( I ) = VAL*TEMPP ( I ) 

10  CONTINUE 
RETURN 
END 
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PSWP.FOR 

PROGRAM  PSWP 

C  THIS  PROGRAM  WILL  READ  THE  OUTPUT  t ILES  OF  XTAL . FOR  AND 
C  GENERATE  A  POWER  SWEEP  AT  A  REQUESTED  FREQUENCY. 

C  IF  THE  FREQUENCY  REQUESTED  IS  ONE  ROW  OF  THE 
C  DATA  FILE  NO  INTERPOLATION  IS  NECESSARY, 

C  IF  THE  FREQUENCY  IS  NOT  THEN  THE  DATA  IS  LINEARLY  INTERPOLATED 
C  BASED  ON  THE  FREQUENCY  DATA  FOR  THE  ROW  ABOVE  AND  BELOW. 

PARAMETER ( I FR= 5 1 4 , I PWR= 1 0 ) 

DIMENSION  FREQ(IFR) , POW ( I FR , I PWR) , VOLT ( IFR , I PWR) 

DIMENSION  VGPH(IPWR) .PGPH(IPWR) .VPLT(IOO) .PPLT(IOO) 

REAL*8  T 1 (IFR) ,T2 (IFR) 

REAL*8  BV(IFR) ,BP(IFR) ,CV(IFR) ,CP(IFR) ,DV(IFR) ,DP(IFR) .VAL1.VAL2 
CHARACTER  LABEL* 40 , FNAME* 20 , XLAB* 10 , YLAB* 10 ,PTITLE*20 
WRITE (*,»)  ’ENTER  THE  VOLTAGE  DATA  FILE  NAME’ 

READ ( * , * )  FNAME 

OPEN ( UNIT= 1 , FILE=FNAME , STATUS = ’OLD* ) 

READd,*)  LABEL 
READ ( 1 , * )  J 
READd,*)  LABEL 

C  GET  THE  NUMBER  OF  ROWS  AND  COLUMNS  IN  THE  DATA  TABLE 
READd,*)  I ROW, ICOL 
READd,*)  LABEL, (Tl( I) ,1*1, ICOL) 

DO  5  1=1 , I ROW 

READd  ,*)  FREQ (I) , (VOLT (I ,J) ,J=1 ,ICOL) 

5  CONTINUE 

10  CLOSE (UNIT= 1 ) 

WRITE (*,*)  ’ENTER  THE  POWER  DATA  FILE  NAME’ 

READ ( * , * )  FNAME 

OPEN (UNIT= 1 , FILE=FNAME , STATUS = ’ OLD ’ ) 

READd,*)  LABEL 
READ (1,«)  J 
READd,*)  LABEL 

C  GET  THE  NUMBER  OF  ROWS  AND  COLUMNS  IN  THE  DATA  TABLE 
READd,*)  IROW.ICOL 
READd,*)  LABEL,  (T1  (I)  ,1  =  1  , ICOL) 

DO  15  1  =  1, 1 ROW 

READ(1 ,  *)  FREQ (I)  ,  (POW (I ,J)  ,J=1  , ICOL) 


READd, »)  LABEL,  (T1  (I)  ,1  =  1  , ICOL) 

DO  15  1=1, I ROW 

READ(1 ,  *)  FREQ (I) , (POW (I ,J) ,J=1 , ICOL) 

CONTINUE 
CLOSE (UNIT=1) 

THE  DATA  TO  BE  ANALYZED  IS  NOW  STORED  IN  THE  ARRAYS  POW  AND  VOLT 

WRITE ( * , * )  ’ENTER  THE  FREQUENCY ( HZ)  THAT  YOU  ARE  INTERESTED  IN’ 
READ ( * , * )  FO 

FIND  OUT  IF  WE  NEED  TO  INTERPOLATE  OR  CAN  WE  USE  THE  DATA  AT  HAND 
IF  (FREQ(l)  . GT . FO . OR . FREQ ( I  ROW)  .LT.FO)  THEN 

WRITE (*,*)  ’THE  FREQUENCY  ’,F0,’  CAN  NOT  BE  USED’ 

STOP 
END  IF 

DO  20  1  =  1  , IROW 

IF  (FREQ (I) . LE . FO . AND . FREQ ( I + 1 ) .GT.FO)  ICH=I 
CONTINUE 

WRITE ( * , * )  FREQ(ICH) , FO , FREQ ( ICH* 1 ) 

IF  (ABS(FREQdCH)  -FO)  .  GT  .  1  .E-2)  THEN 


LE . FO . AND . FREQ ( I + 1 ) .GT.FO)  ICH=I 


v  v/V  vW. 
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C  TIME  TO  INTERPOLATE 
C  INTERPOLATE  THE  VOLTAGE 
DO  25  1*1  ,  ICOT. 

T I (1) =FREQ ( ICH) 

T 1 (2) =  FREQ ( ICH+ 1 ) 

T2 ( 1 ) = VOLT ( ICH , I ) 

T2 (2) =VOLT ( ICH+ 1 ,1) 

N=2 

CALL  SPCOEF (N , T 1 , T2 , BV , CV , DV , I FLAG) 

IF  (IFLAG.LT. 0)  THEN 

WRITE ( * , * )  ’ERROR  IN  VOLTAGE  INTERPOLATION’ 
WRITE ( * , # )  IFLAG 
STOP 
ENDIF 
VAL 1 =F0 

CALL  SVALUE (N . T 1 , T2 , BV , CV , DV , VAL1 , VAL2 , IFLAG) 
IF  (IFLAG.LT. 0)  THEN 

WRITE ( * , # )  ’ERROR  IN  EVALUATION’ 

WRITE (  *  ,  * )  IFLAG 
STOP 
ENDIF 

VGPH(I) =VAL2 
25  CONTINUE 

C  INTERPOLATE  THE  POWER 
DO  30  1=1 ,ICOL 

T 1 ( 1 ) =FREQ ( ICH) 

T 1 (2) =  FREQ ( ICH+ 1 ) 

T2 ( 1 ) *POW ( ICH , I ) 

T2 (2) =P0W(ICH+1 , I) 

N=2 

CALL  SPCOEF (N , T 1 , T2 , BP , CP , DP , IFLAG) 

IF  ( IFLAG. LT.O)  THEN 

WRITE  (.*,*)  'ERROR  IN  POWER  INTERPOLATION’ 
WRITE  {.*,*)  IFLAG 
STOP 
ENDIF 
VAL 1 =F0 

CALL  SVALUE (N,T1 , T2 , BP , CP , DP , VAL 1 ,VAL2, IFLAG) 
IF  (IFLAG. LT.O)  THEN 

WRITE ( « , » )  ’ERROR  IN  POWER  EVALUATION’ 

WRITE ( * , * )  IFLAG 
STOP 
ENDIF 

PGPH ( I ) = VAL2 
30  CONTINUE 

ELSE 

DO  35  1=1 . ICOL 

VGPH(I) =VOLT ( ICH , I ) 

PGPH ( I ) =POW( ICH , I ) 

35  CONTINUE 

ENDIF 
C 

C  LET’S  INTERPOLATE  ON  PGPH  AND  VGPH  TO  GET  A  SMOOTH  CURVE 
DO  40  J= 1 , I PWR 

T  1  ( J) =DBLE (VGPH(J) ) 
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T2(J) =DBLE ( PGPH ( J) ) 

CONTINUE 

CALL  SPCOEF ( ICOL , T 1 , T2 , BP , CP , DP , I FLAG) 


STEP=  (VGPH(ICOL)  -,,GPH{  1)  )  / 99 
DO  45  J=  1  , 100 

VPLT(J) = VGPH ( 1 ) + (J- 1) *STEP 
VAL1=DBLE (VPLT ( J) ) 

CALL  SVALUE ( ICOL , T 1 , T2 , BP , CP , DP , VAL 1 , VAL2 . IFLAG) 
PPLT(J) =SNGL ( VAL2 ) 

45  CONTINUE 

C  LET’S  HAVE  THE  COMPUTER  EVALUATE  SOME  XTAL  VOLTAGES  INTO 
C  BASED  ON  THE  INTERPOLATED  DATA  POINTS 

50  WRITE ( * , * )  'ENTER  THE  XTAL  VOLTAGE (VOLTS) ' 

READ (  *  ,  * )  2VAL 
VAL 1 =DBLE (2VAL) 

CALL  SVALUE (ICOL ,T1 , T2 , BP , CP , DP , VAL 1 , VAL2 , IFLAG) 

WRITE ( » ,51)  ZVAL , SNGL ( VAL2 ) .IFLAG 

51  FORMAT ( ’  ’ , ’ FOR  ’.1PE10.3,’  VOLTS,  ’.1PE10.3,’  WATTS’, 14) 

52  WRITE («,»)  ’FOR  MORE  DATA  ENTER  A  1’ 

WRITE («,«)  ’TO  CONTINUE  ENTER  A  0’ 

WRITE («,«)  ’TO  STOP  ENTER  A  2’ 

READ ( * , * )  ICNTR 
IF  ( ICNTR. EQ.l)  GOTO  50 
IF  (ICNTR. EQ. 2)  STOP 
IF  (ICNTER.NE.O)  GOTO  52 

C  WE  ARE  NOW  READY  TO  PLOT  PGPH  VS  VGPH 
C  USE  VAX  GRAPHICS 

C  GRAPHICS  FOR  USE  ON  4010,4014,4025,4105 
WRITE (*,*)  ’ENTER  THE  DEVICE  ID* ’ 

READ ( » , » )  IDEVIC 

IF  (IDEVIC. EQ. 4014. OR. IDEVIC. EQ. 4105)  THEN 
XMI =.01 
XMA= .9 
YMI =.01 
YMA= .9 

ELSEIF  (IDEVIC. EQ. 4010)  THEN 
XMI =.01 


POWERS 


(IDEVIC. EQ. 4025) 
.  05 


XMI = . 0 1 
XMA= .9 
YMI =.01 
YMA= .9 

ELSEIF  (IDEVIC. EQ. 4010)  THEN 
XMI =.01 
XMA= . 6 
YMI = . 0 1 
YMA= . 6 

ELSEIF  (IDEVIC. EQ. 4025)  THEN 
XMI =.05 
XMA= .85 
YMI = .05 
YMA= .85 
IDEVIC=40 10 
ELSE 

WRITE ( * , » )  ’DEVICE’ , IDEVIC , ’  NOT  USED’ 
STOP 
END  IF 

WRITE (*.»)  ’ENTER  A  1  FOR  A  LINEAR  PLOT’ 
WRITE ( » , » )  ’ENTER  A  2  FOR  A  LOG-LOG  PLOT’ 
READ ( * , * )  IOPT 

IF  (IOPT.NE.l  .AND.  I0PT.NE.2)  GOTO  100 


S  vV  A  a  .k  .V  A- .%  .v,  •  A  A 
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CALL  GSTART ( IDEVIC) 

CALL  NEWPAG 

CALL  MINMAX ( VGPH , I COL , XWMIN , XWMAX) 

CALL  MINMAX (PGPH , ICOL , YWMIN , YWMAX) 

IF  (IOPT.EQ.l)  THEN 

CALL  MAP (XWMIN, XWMAX, YWMIN, YWMAX, XMI , XMA , YMI , YMA) 
CALL  XLINCXWMIN, YWMIN, XWMIN, XWMAX, 5) 

CALL  YLIN (XWMIN, YWMIN, YWMIN, YWMAX, 5) 

ELSEIF  (IOPT.EO.2)  THEN 

CALL  MAPLL( XWMIN, XWMAX, YWMIN, YWMAX, XMI , XMA, YMI , YMA) 
CALL  XLOG( XWMIN, YWMIN, XWMIN, XWMAX) 

CALL  YLOG ( XWMI N , YWMI N , YWMI N , YWMAX) 

ENDIF 

XLAB= 'VOLTS’ 

YLAB= ’POW(W) ’ 

PTITLE= ’ POWER  SWEEP’ 

CALL  XLABEL (XLAB) 

CALL  YLABEL ( YLAB ) 

CALL  TITLE (PTITLE) 

CALL  CURVE (VGPH, PGPH. ICOL, 0,-1) 

CALL  CURVE (VPLT.PPLT, 100,0,0) 

CALL  GRSTOP 


